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ABSTRACT 


A  7-ft  cubical  sample  of  dry  sand  was  tested  using  the  tri- 
axial  testing  device  constructed  by  Kopperman  et  al  (1982)  and  Knox  et  al 
(1982).  The  sand  was  the  same  as  that  used  by  Kopperman  and  Knox.  A 
new  raining  device  was  fabricated  and  used  to  construct  this  sand  sample 
which  resulted  in  a  more  uniform  sample  than  prepared  earlier.  Improve¬ 
ments  were  also  made  to  the  excitation  ports  in  order  to  have  better 
control  of  the  stress  states. 

Extensive  tests  were  performed  under  the  following  different 
stress  states:  isotropic,  biaxial  and  triaxial.  In  each  case,  velo¬ 
cities  of  P-waves  propagating  along  all  principal  stress  directions  were 
measured.  Results  from  these  tests  lead  to  the  following  conclusions: 

1.  the  effect  of  stress  history  on  P-wave  velocity  is  negligible,  2.  the 
sample  can  be  treated  as  a  cross-anisotropic  material  under  isotropic 
confinement  due  to  structural  anisotropy,  3.  complete  anisotropy 
resulted  by  the  coupling  of  stress  anisotropy  and  structural  anisotropy, 
and  4.  P-wave  velocity  depends  on  the  principal  effective  stress  in  the 
direction  of  propagation,  with  principal  stresses  perpendicular  to  the 
direction  of  propagation  having  a  negligible  effect  on  velocity. 
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CHAPTER  ONE 


INTRODUCTION 


A  large-scale  trlaxlal  testing  device  was  designed  and  con¬ 
structed  during  1980  and  1981  under  the  sponsorship  of  a  grand  from  the 
United  States  Air  Force  Office  of  Scientific  Research  (AFOSR).  With 
this  device,  laboratory  studies  of  the  effects  of  Isotropic,  biaxial 
and  trlaxlal  states  of  stress  on  the  velocity  of  compression  and  shear 
waves  propagating  through  soil  can  be  performed.  An  initial  series  of 
tests  with  dry  sand  was  performed  by  Kopperman  et  al  (1982)  on  com¬ 
pression  waves  and  by  Knox  et  al  (1982)  on  shear  waves.  The  objectives 
of  the  research  reported  herein  are:  1.  to  check  the  validity  of  the 
conclusions  reported  by  Kopperman  et  al  on  compression  waves 
(P-waves),  and  2.  to  extend  this  earlier  work  by  performing  additional 
tests  under  anisotropic  states  of  stress. 

1.1  BACKGROUND  INFORMATION 

The  trlaxlal  testing  device  is  essentially  a  reinforced 
steel  box  with  Interior  dimensions  of  7  ft  (2.1  m)  on  each  side. 
Sketches  of  the  device  are  shown  In  Figs.  1.1  and  1.2,  and  the  com¬ 
pleted  device  Is  shown  In  Fig.  1.3.  Equipment  associated  with  the 
device  Is  used  to  :  1.  place  sand  Into  the  device,  2.  pressurize  the 


Fig.  1.3  -  Completed  Triaxial  Device  (From 
Kopperman  et  al,  1982) 


sand  mass  to  the  desired  stress  state,  3.  generate  compression  or  shear 
waves  In  the  sand  mass,  4.  monitor  and  digitally  record  these  wave¬ 
forms,  and  5.  monitor  stress  and  strain  throughout  the  sand  sample  dur¬ 
ing  testing.  A  schematic  drawing  of  the  trlaxlal  device  and  associated 
systems  Is  shown  In  Fig.  1.4. 

A  new  raining  device  was  designed  and  constructed  as  part  of 
the  work  reported  herein  (Section  2.3).  This  was  done  in  an  attempt  to 
place  the  sand  Into  the  trlaxlal  device  In  a  more  uniform  state  than 
possible  by  Knox  and  Kopperman.  With  the  new  rainer,  the  rate  of  sand 
flowing  into  the  device  during  sample  construction  was  carefully  con¬ 
trolled.  As  a  result,  the  density  of  the  sand  specimen  increased  by 
about  six  percent  compared  to  the  earlier  tests  (Section  3.3),  and  the 
specimen  was  very  uniform  with  a  maximum  variation  in  density  of  less 
than  6.3  percent. 

Each  wall  of  the  trlaxlal  device  Is  designed  to  represent  a 
principal  plane  so  that  axes  perpendicular  to  the  walls  of  the  device 
represent  principal  directions.  To  permit  Independent  control  of  the 
pressure  In  each  of  the  three  principal  directions,  confining  pres¬ 
sures  are  applied  to  the  sand  mass  using  three  membranes  (water  pres¬ 
sure  bags)  placed  on  the  Inside  of  the  trlaxlal  device  :  one  on  the  top 
and  one  on  each  of  two  adjacent  sides,  namely  the  north  and  west  sides 


of  the  device.  Each  membrane  has  two  ports  located  at  opposite  corners 
so  that  water  can  be  allowed  to  fill  up  or  drain  from  the  membrane.  One 
of  these  membrane  Is  shown  In  Fig.  1.5.  When  the  membranes  full  of 
water,  air  pressure  from  the  building  air  supply  Is  used  to  pressurize 
the  membranes.  The  air  pressure  Is  monitored  using  air  regulators 
together  with  a  12-1 n.  (30.5  cm)  diameter,  Helse  type  CM  pressure  gauge 
(accurate  to  within  ±0.1  percent  of  full-scale  reading).  A  schematic 
drawing  of  the  pressurizing  system  Is  shown  In  Fig.  1.6,  and  a  picture 
of  the  control  panel  Is  shown  In  Fig.  1.7. 

Since  the  overall  objective  of  this  research  Is  to  study  body 
waves  (compression  and  shear  waves)  propagating  through  a  soil  mass,  a 
mechanism  for  generating  these  waves  at  the  soil  boundaries  Is  neces¬ 
sary.  To  achieve  this  goal,  excitation  ports  were  placed  in  the  center 
of  the  three  mutually  perpendicular  sides  of  the  trlaxlal  device  which 
were  not  covered  by  membranes  :  the  bottom,  and  the  south  and  east 
sides.  The  excitation  ports  were  modified  and  Improved  from  those  used 
by  Kopperman  et  al  (1982)  as  discussed  In  Section  2.4.  (A  detailed 
sketch  of  an  excitation  port  Is  shown  In  Fig.  2.5.) 

Compression  and  shear  waves  propagating  through  the  sand 


mass  are  monitored  and  recorded  with  the  electronics  shown  schemat- 


(From  Kopperman  et  al,  1982) 


On/Off  Toggle  Valve 


(From  Kopperman  et  al,  1982) 


leal ly  In  Fig.  1.8.  Improvements  and  recalibration  of  Instruments  In 
this  monitoring  and  recording  system  are  discussed  In  Section  2.5. 
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1.2  PRESENT  TESTING 

Once  all  modifications  and  calibrations  were  completed, 
tests  with  various  states  of  confinement  were  performed.  The  proce¬ 
dures  and  results  for  Isotropic,  biaxial  and  trlaxlal  confinement 
states  are  presented  In  Chapters  4,  6  and  7,  respectively. 

The  effects  of  stress  history  and  structural  anisotropy  on 
P-wave  velocity  were  determined  first  before  any  conclusion  about  the 
effects  of  confinement  on  P-wave  velocity  could  be  reached.  These 
effects  are  most  easily  seen  In  the  results  presented  In  Chapter  4  for 
Isotropic  confinement.  The  effect  of  stress  history  on  P-wave  velocity 
Is  small  and  can  be  neglected.  Due  to  structural  anisotropy,  the  sand 
sample  under  Isotropic  confinement  can  be  treated  as  a 
cross-anisotropic  material  with  P-wave  velocities  In  the  horizontal 
directions  equal  and  about  17  percent  higher  than  the  velocity  In  the 
vertical  direction.  A  mathematical  treatment  for  the  stress-strain 
relationship  In  cross-anisotropic  model  Is  presented  In  Chapter  5. 


Monitoring  and  Recording  Systems 


The  results  of  P-wave  velocity  measured  under  Isotropic, 
biaxial  and  trlaxlal  confinement  states  lead  to  the  same  conclusion: 
that  Is,  P-wave  velocity  Is  proportional  to  the  effective  principal 
stress  In  the  direction  of  wave  propagation  raised  to  a  positive  power 
less  than  unity.  The  normal  principal  stresses  perpendicular  to  the 
direction  of  wave  propagation  have  a  negligible  effect  on  P-wave  veloc¬ 
ity.  This  confirms  the  earlier  work  by  Kopperman  et  al  (1982)  In  a  more 
complete  fashion. 


CHAPTER  TWO 

IMPROVEMENTS  OF  TESTING  FACILITIES 

2.1  INTRODUCTION 

Considerable  effort  was  expended  on  this  project  to  Improve 
the  trlaxlal  loading  device  and  associated  systems  used  by  Kopperman  et 
al  (1982)  and  Knox  et  al  (1982).  This  work  was  done  so  that  duplicate 
tests  could  be  performed  under  more  carefully  controlled  conditions 
than  possible  In  the  earlier  tests.  In  addition,  larger  variations  In 
tests  under  anisotropic  loading  were  performed  than  done  earlier 
because  of  time  limitations.  Improvements  were  made  In  the  following 
areas  :  1.  construction  of  a  storage  facility  In  which  the  sand  could 
be  kept  dry,  2.  development  of  a  new  raining  device  used  to  build  spec¬ 
imens,  3.  addition  of  the  capability  to  control  the  load  applied  at 
each  excitation  port,  4.  calibration  of  the  accelerometers,  stress 
cells  and  strain  cells,  and  5.  Inclusion  of  additional  accelerometers 
to  provide  more  complete  data  for  an  attenuation  study  and  for  an  Ini¬ 
tial  study  of  shear  waves  which  are  polarized  obliquely  (at  some  angles 
other  than  parallel  to  the  principal  stress  axes). 


2.2  STORAGE  OF  SAND 


One  of  the  time-consuming  problems  encountered  In  the  earli¬ 
er  study  was  drying  of  moist  sand  (which  was  stored  outdoors)  before 
sample  construction.  Therefore,  It  was  decided  to  construct  some  type 
of  storage  facility.  A  corrugated-steel  bln  measuring  about  six  ft 
(1.83  m)  In  diameter  and  six  ft  (1.83  m)  In  height  was  available  at  no 
cost  in  Cockrell  Hall.  The  bln  was,  therefore,  moved  to  the  Bal cones 
Research  Center  and  was  located  outside  the  Phil  M.  Ferguson  Struc¬ 
tural  Laboratory,  next  to  the  west  wall  and  close  to  the  north  truck 
entrance  of  the  Laboratory.  This  Is  a  convenient  location  since  the 
trlaxlal  device  Is  located  at  the  north-west  corner  Inside  the  Labora¬ 
tory. 


To  keep  the  sand  In  the  bln  dry,  the  bln  was  placed  on  top  of 
a  steel  table  as  shown  In  Fig.  2.1.  The  top  of  the  table  was  made  of  a 
6.33-ft  (1.93  m)  square  plate  with  a  thickness  of  0.25  In.  (0.64  cm). 
Four  legs  and  a  framework  for  the  table  top  were  built  of  WI8X10 
I-beams.  The  legs  were  2-ft  (61.0  cm)  high  and  were  welded  to  the 
framework  at  the  top  and  to  0.5-1n.  (1.27  cm)  thick,  1-ft  (30.5  cm) 
square  steel  plates  at  the  bottom.  The  bottom  plates  were  used  to  dis¬ 
tribute  load  from  each  leg  to  a  concrete  footing  which  was  also  1-ft 
(30.5  cm)  square  and  about  1.5-ft  (45.7  cm)  deep. 


The  lower  edge  of  the  bln  was  fixed  to  the  top  of  the  steel 
table  with  small  pieces  of  angle  iron  with  bolt  connections.  An  open¬ 
ing  of  about  8  in. (20. 3c  cm)  by  11  in.  (27.9  cm)  was  cut  near  the  bottom 
edge  of  the  bin.  A  trap  door  with  hinges  at  the  top  was  bolted  to  the 
bin  to  cover  the  opening.  Styrofoam  was  glued  under  the  perimeter  ->f 
the  door  to  prevent  sand  leakage  and  to  provide  a  moisture  seal.  The 
trap  door  could  be  opened  to  any  desired  position  with  the  control  of  a 
bolt  connection.  Once  the  trap  door  was  opened,  dry  sand  would  flow  out 
of  the  bin  into  a  concrete  bucket  seated  next  to  it.  This  greatly 
improved  the  sample  construction  process  (Section  3.3). 

A  6.5-ft  (1.9  m)  square  top  was  built  for  the  bln  with  ply¬ 
wood  covered  by  sheet  metal.  It  was  decided  that  this  was  the  most  eco¬ 
nomical  way  to  provide  a  waterproof  top  for  the  bln.  Two  handles  were 
bolted  to  the  top  so  that  It  could  be  removed  and  replaced  easily  with 
the  use  of  a  fork-lift. 

2.3  SOIL  PLACEMENT 

Since  the  old  ralner  constructed  by  Knox  and  Kopperman  could 
not  control  the  amount  the  hopper  gate  opened  and  hence  could  not  con¬ 
trol  the  flow  rate  of  the  sand  (Kopperman  et  al ,  1982),  a  new  ralner 


was  designed  and  constructed.  The  ralner  is  simply  an  open  rectangular 
box  measuring  7.25-ft  (2.21  m)  long,  1.75-ft  (53.3  cm)  wide  and  1.5-ft 
(45.7  cm)  deep.  The  box  Is  made  of  plywood  and  has  a  steel  framework  on 
the  outside  as  shown  In  Figs.  2.2  and  2.3.  The  plywood  bottom  has  four 
rows  of  0.75-in.  (1.91  cm)  diameter  holes  drilled  In  It  In  the  pattern 
shown  In  Fig.  2.4.  Each  end  of  the  ralner  has  a  3.42-ft  (104.1  cm)  long 
angle  (L3X3)  welded  to  the  lower  edge  of  the  steel  framework.  Heav¬ 
y-duty  castors  are  fixed  to  each  end  of  the  L3X3  angle.  These  extended 
portions  with  castors  add  to  the  stability  of  the  ralner  when  It  Is 
being  moved  during  the  raining  operation.  Four  trap  doors  were  built 
under  the  plywood  bottom,  each  covering  two  rows  of  holes  as  shown  in 
Fig.  2.4.  The  amount  the  trap  door  is  opened  Is  controlled  by  a 
lever-arm  system.  The  lever-arm  system  can  be  securely  fixed  with  the 
trap  doors  closed  or  with  the  trap  doors  opened  to  the  selected  posi¬ 
tion  to  allow  the  desired  flow  rate  of  sand.  A  wire  mesh  with  openings 
of  0.25  In.  (6.4  mm)  was  placed  below  the  ralner  as  shown  In  Fig.  2.3. 
The  wire  mesh  acts  as  a  dispersing  screen  since  sand  striking  the  mesh 
will  bounce  off  randomly,  hence  constructing  a  more  uniform  specimen. 
The  screen  also  keeps  foreign  materials  like  gravels  and  organics  from 
becoming  part  of  the  sand  sample. 

The  same  wooden  collar  used  by  Knox  and  Kopperman  was  used  to 
obtain  the  same  range  of  drop  height  of  the  sand  Into  the  trlaxlal 


Legend:  1.  Steel  framework  of  the  new  rainer 

2.  Rail  along  which  rainer  travels 

3.  Angle  Irons  (L2X2)  along  the  corners  of  the  wooden  collar 

4.  Wooden  collar 

5.  Castor  of  rainer 


Fig.  2-2  View  of  New  Rainer  Loaded  with  Sand  and  in  Place  for  Raining 
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Legend:  1.  New  rainer 

2.  Control  Arm  (see  Fig.  2.4) 

3.  Lever  connection  to  trap  doors 

4.  Wire  mesh  used  as  dispersing  screen 

Fig.  2.3  Raining  Sand  into  Trlaxial  Device  Using  New  Rainer 


New  Rainer 
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device,  namely  9.5  ft  (2.9  m)  at  the  start  of  the  filling  operation  to 
2.5  ft  (76.2  cm)  at  the  conclusion.  Since  the  new  ralner  is  much  heavi¬ 
er  than  the  old  one,  the  collar  was  reinforced  with  angle  irons  (L2X2) 
along  the  four  vertical  edges.  Angle  irons  at  two  corners  are  shown  in 
Fig.  2.2.  The  angles  were  bolted  to  the  collar  and  are  3-ft  (91.4  cm) 
in  length  so  that  they  support  the  rails  (as  shown  in  Fig.  2.2)  at  the 
top  while  resting  on  the  top  longitudinal'  reinforcement  ring  of  the 
triaxial  device.  This  greatly  increased  the  rigidity  and  stability  of 
the  collar.  Each  rail  also  has  three  sections  of  6-in.  (15.2  cm)  long 
angles  (L2X2)  bolted  to  the  collar  and  distributed  evenly  under  the 
rail  to  provide  better  support  for  the  rail. 

2.4  EXCITATION  PORTS 


Sand  leakage  from  the  two  horizontal  excitation  ports  (see 
Fig.  1.3)  was  reported  In  the  previous  tests.  This  was  most  probably 
due  to  a  slight  rotation  of  the  anvil  and  subsequent  sliding  of  sand 
out  of  a  gap  formed  between  the  anvil  and  hole  In  the  wall  for  the  exci¬ 
tation  port.  To  prevent  this,  a  clamping  tool  was  designed  to  hold  the 
anvil  in  place  when  testing  was  not  being  performed.  The  clamping  tool 
Is  basically  two  similar  pieces  of  steel  bars  placed  together.  Three 
holes  have  been  drilled  in  the  bars  so  that  the  shaft  of  the  anvil  fits 
in  the  center  hole  and  the  two  large  bolts  of  the  external  frame  fit  in 
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each  of  the  side  holes.  The  clamping  tool  Is  shown  In  place  on  an  exci¬ 
tation  port  In  Fig.  2.5.  The  ends  of  the  two  steel  bars  can  be  bolted 
tightly  together,  hence  the  anvil  Is  held  fixed  with  respect  to  the 
external  frame,  which  In  turn  Is  fixed  to  the  trlaxlal  device.  The 
anvil  In  each  excitation  port  was  clamped  with  these  tools  except  dur¬ 
ing  testing. 

It  was  decided  that  the  anvil  of  the  shaft  would  provide 
enough  striking  surface  to  generate  the  desired  waves,  therefore  the 
3-in.  (7.6  cm)  square  plate  used  previously  as  the  striking  surface  was 
removed.  Shear  waves  were  generated  in  the  soil  by  striking  the  anvil 
shaft  either  horizontally  or  vertically  (parallel  to  the  side  of  the 
trlaxlal  device)  about  1  In.  (2.5  cm)  from  the  wall.  Compression  waves 
were  generated  by  striking  the  top  surface  of  the  anvil  in  the  direc¬ 
tion  of  the  axis  of  anvil  (perpendicular  to  the  side  of  the  trlaxial 
device).  The  striking  locations  are  shown  in  Fig.  2.5b. 

The  anvil  In  each  excitation  port  which  Is  In  contact  with 
the  soil  must  maintain  the  same  pressure  against  the  soil  as  Is  being 
applied  by  the  membrane  on  the  opposite  face  of  the  trlaxlal  device.  To 
monitor  the  load  applied  by  the  anvil,  strain  gages  were  mounted  to  the 
thin  plate  of  the  external  frame  of  the  port  as  shown  In  Fig.  2.5b. 
After  a  preliminary  calculation  on  the  range  of  bending  strain  of  the 
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thin  plate  under  the  working  stress  range  of  10  to  40  psi  (68.9  to 
275.6  kPa),  strain  gages  of  type  EA-06-500BL-350  were  selected.  This 
is  a  foil-type,  bonded,  resistive  gage  with  350  ohms  ±0.15  percent  as 
the  gage  resistance  and  2.07  ±0.5  percent  as  the  gage  factor.  The 
strain  gages  and  Installation  aids  were  purchased  from  the 
Micro-Measurements  Division  of  the  Measurements  Group  at  Raleigh, 
North  Carolina. 

The  strain  gages  were  mounted  on  the  thin  plate  following  the 
procedures  recommended  by  the  manufacturer.  The  relation  between 
strain  and  resistance  change  for  these  gages  Is  : 

£=(1/F)*(AR/R) 

where:  e  =  strain  In  microstrains  (10~*  cm/cm), 

F  =  gage  factor, 

AR  =  change  In  resistance  In  ohms,  and 
R  =  gage  resistance  In  ohms. 

A  strain  indicator  was  used  to  read  off  the  resulting  strain 
digitally  in  microstrains.  Two  strain  gages  were  mounted  to  each  plate 
(one  on  top  and  one  beneath)  In  such  a  pattern  that  temperature  compen¬ 
sation  was  provided  and  the  output  was  doubled.  After  the  gages  were 
connected  to  the  lead  wires,  silicone  glue  was  placed  over  the  gages  to 
protect  them  from  the  ambient  environment  of  moisture,  dust  and  dirt. 
A  completed  excitation  port  is  shown  In  Fig.  2.6. 


The  strain  gages  used  as  secondary  transducers  were  cali¬ 
brated  In  terms  of  pressure  on  the  face  of  the  anvil  In  contact  with 
the  sand.  Therefore,  the  strain  gages  were  calibrated  using  a  load  cell 
to  monitor  the  applied  force.  The  area  of  the  anvil  In  contact  with  the 
sand  Is  9  sq-ln.  (522.6  cm2),  hence  a  1-lb  force  (4.45  N)  Induces  a 
pressure  of  0.111  psl  (0.766  kPa)  under  the  face  of  the  anvil.  The 
calibration  set-up  Is  shown  In  Fig.  2.7.  A  loading  frame  was  used  to 
push  the  anvil  downwards  (analogous  to  the  soil  pushing  the  anvil  out¬ 
wards)  and  a  load  cell  was  seated  between  the  anvil  and  the  loading 
piston  of  the  press  to  monitor  the  applied  pressure.  Unfortunately  at 
the  time  of  calibration,  only  a  300-lb  (1.34  kN)  load  cell  was  avail¬ 
able.  Therefore,  the  maximum  reading  corresponded  to  a  pressure  of  33.3 
psl  (229  kPa),  which  Is  slightly  less  than  the  maximum  pressure  of  40 
psl  (275.6  kPa)  used  In  these  tests.  Two  strain  Indicators  were  used: 
one  to  monitor  the  load  cell  and  the  other  the  strain  gages.  The  final 
calibration  curves  are  shown  In  Figs.  2.8  to  Fig.  2.10.  The  reading 
for  40  psl  (275.6  kPa)  was  obtained  by  linear  extrapolation  since  the 
material  Is  essentially  linearly  elastic. 

The  calibration  curves  In  Figs.  2.8  to  2.10  showed  small  hys¬ 
teresis  upon  load  cycling.  To  see  If  creep  had  an  Insignificant  effect 
on  the  strain  gages,  the  top-bottom  (TB)  excitation  port  was  kept  under 


Legend:  1.  300-lb  load  cell 

2.  Strain  gages 

3.  Strain  Indicator  for  strain  gages 

4.  Strain  Indicator  for  load  cell 

Fig.  2.7  Set-Up  for  Calibrating  Strain  Gages  on  Each  Excitation  Port 


Fig.  2.9  Calibration  Curve  for  Strain  Gages  on 
NS  Excitation  Port 


Land  Cell  Reading  (lbs) 


Fig.  2.10  Calibration  Curve  for  Strain  Gages  on 
EW  Excitation  Port 


a  pressure  of  13.8  psl  (95.1  kPa)  for  a  period  of  ten  days.  Then  the 
strain  gages  were  calibrated  again.  The  data  essentially  fell  within 
the  hysteresis  loops  as  shown  In  Fig.  2.8,  and  therefore  creep  was 
assumed  to  be  negligible. 

As  an  example  of  the  use  of  the  calibration  curves,  assume  an 
Isotropic  pressure  of  15  psl  (103.4  kPa)  was  applied  to  the  sand  speci¬ 
men.  The  corresponding  average  strain  reading  on  the  TB  port  should 
then  be  152  microstrains  as  shown  In  Fig.  2.8.  In  order  for  the  anvil 
contacting  the  sand  to  attain  this  pressure,  the  adjustment  screw  (Fig. 
2.5)  In  the  excitation  port  would  be  either  tightened  or  loosened  until 
the  strain  reading  of  the  strain  gages  was  152  microstrains. 

2.5  MONITORING  AND  RECORDING  SYSTEM 


The  same  Instruments  (accelerometers,  charge  amplifiers, 
oscilloscopes,  stress  cells  and  strain  cells)  used  In  the  earlier  tests 
were  used  In  the  monitoring  and  recording  system,  except  that  a  signif¬ 
icant  effort  was  made  to  recalibrate  the  instruments  before  this  series 
of  dynamic  tests  was  performed.  In  addition,  eight  more  accelerometers 
were  placed  Inside  the  trlaxlal  device  making  a  total  of  29  accelerome¬ 
ters.  All  accelerometers  were  placed  along  the  principal  axes  of  the 
sand  mass.  Four  of  the  eight  new  accelerometers  were  added  along  the 


vertical  direction  to  monitor  vertically  propagating  shear  waves  which 
were  polarized  obliquely,  and  the  other  four  accelerometers  were  added 
along  the  north-south  direction  to  provide  more  data  for  the  atten¬ 
uation  study. 

2.5.1  Digital  Oscilloscopes 

Two  digital  oscilloscopes  with  magnetic  storage  capabili¬ 
ties,  series  2090,  were  purchased  earlier  from  the  Nlcolet  Instrument 
Corporation  at  Madison,  Wisconsin.  These  units  were  used  to  monitor 
and  record  the  outputs  from  the  accelerometers.  One  of  the  two  oscil¬ 
loscopes  (Serial  No.  79703)  was  sent  back  to  the  dealer  for  mainte¬ 
nance,  cleaning  and  adjustment  of  the  magnetic  storage  unit  before  use. 

2.5.2  Accelerometers 

Twenty-one  Endevco  Isoshear  accelerometers,  model  7701-100, 
were  used  In  the  previous  tests.  These  accelerometers  have  a  typical 
charge  sensitivity  of  100  pc/g  (10~12  coulombs  per  gravitational  accel¬ 
eration).  Two  more  accelerometers  of  this  model  and  six  accelerometers 
of  model  7701-50  were  added  In  the  present  tests.  All  model  7701-50 
accelerometers  have  half  the  charge  sensitivity  of  a  model  7701-100, 
namely  50  pc/g.  This  variation  In  charge  sensitivity  did  not  affect 


collection  of  data  since  the  same  full-scale  output  could  be  obtained 
simply  by  changing  the  sensitivity  dial  and  full-scale  range  switch  on 
the  charge  amplifier,  which  Is  used  to  amplify  the  signal  before  it  is 
displayed  on  the  oscilloscopes.  The  six,  7701-50,  accelerometers  were 
used  In  the  locations  of  the  first  and  last  3-D  accelerometer  packages 
In  the  east-west  direction  as  shown  in  Fig.  2.11. 

Since  the  accelerometers  are  precise  instruments  which  are 
designed  to  give  many  years  of  service  without  maintenance,  only  a  com¬ 
parative  technique  was  used  to  calibrate  the  accelerometers.  Any  one 
accelerometer  was  believed  to  be  reliable  and  could  be  picked  as  a 
"standard"  against  which  the  other  accelerometers  could  be  compared. 
As  such,  the  calibration  procedure  consisted  of  Imposing  "identical" 
motions  on  the  subject  and  standard  accelerometers  and  the  two  outputs 
were  compared.  It  was  found  that  all  accelerometers  were  functioning 
satisfactory  with  differences  between  outputs  from  the  subject  and 
standard  accelerometers  less  than  ±0.05  percent. 

2.5.3  Accelerometer  Blocks 

Four  new  accelerometer  blocks  were  added  In  the  trlaxlal 
device.  Each  new  accelerometer  block  was  designed  to  house  two  accel¬ 
erometers  as  compared  to  the  3-D  accelerometer  blocks  constructed  ear- 
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Legend:  1.  NS  =  north-south  principal  axis 
EW  =  east-west  principal  axis 
TB  *  vertical  (top-bottom)  principal  axis 

2.  Letter  following  accelerometer  number  indicates  direction 
of  sensitivity  for  that  particular  accelerometer 
(v*vertical,  s=south  and  w=west) 

3.  Accelerometers  23,  25  are  45°  to  N-S  axis 
Accelerometers  24,  26  are  22.5°  to  N-S  axis 

4.  SN-7,  SN-8  are  two  pairs  of  4"  strain  sensors  for 
vertical  strain  measurement 

Fig.  2.11  -  Schematic  View  of  Instrumentation  Locations 


Her  which  housed  three  accelerometers.  A  schematic  drawing  of  the 
location  of  the  accelerometer  blocks  Inside  the  trlaxial  device  is 
shown  In  Fig.  2.11.  Two  new  accelerometer  blocks  were  located  In  the 
north-south  (NS)  principal  direction  so  that  more  data  could  be  col¬ 
lected  In  that  direction  for  an  attenuation  study  in  another  phase  of 
the  project  being  conducted  by  Dr.  Carroll  at  West  Point.  The  two  other 
new  accelerometer  blocks  were  located  in  the  vertical  principal  axis  to 
monitor  vertically  propagating  shear  waves  which  were  polarized 
obi Iquely. 

One  of  the  new  accelerometer  blocks  was  machined  from  alumi¬ 
num  while  the  other  three  were  made  of  Birch  wood,  (the  same  material 
as  that  of  the  3-D  blocks).  A  picture  of  the  aluminum  block  Is  shown  In 
Fig.  2.12.  The  aluminum  block  was  attached  to  the  excitation  port  on 
the  south  wall  of  the  triaxial  device.  The  area  of  the  square  face  of 
the  block  Is  same  as  the  base  of  anvil  to  which  the  block  is  rigidly 
bolted.  One  accelerometer  inside  the  block  Is  oriented  In  the 
north-south  direction  and  the  other  east-west  (accelerometers  Is  and 
2w,  respectively.  In  Fig.  2.11).  These  accelerometers  were  used  to  pro¬ 
vide  more  data  for  the  attenuation  study  being  conducted  by  Dr.  Car- 
roll.  Two  more  accelerometers  oriented  In  the  same  relative  directions 
(accelerometers  6s  and  7w  In  Fig.  2.11)  and  rigidly  attached  to  a  new 
wooden  block  were  located  2.5-ft  (1.07  m)  away  from  the  aluminum  block 


Fig.  2.12  Accelerometer  Being  Assembled  in  Aluminum  Accelerometer 
Block  Which  is  Part  of  NS  Excitation  Port 


along  the  north-south  principal  direction.  These  accelerometers  were 
also  Included  for  use  In  the  attenuation  study.  These  new  wooden  blocks 
are  of  the  same  size  as  the  3-D  blocks. 

The  other  new  wooden  blocks  were  added  in  the  vertical  direc¬ 
tion:  one  was  located  1-ft  (30.5  cm)  above  the  central  accelerometer 
block  and  the  other  1-ft  (30.5  cm)  below  as  shown  in  Fig.  2.11.  Each 
new  wooden  block  has  holes  for  housing  the  accelerometers  machined  at 
different  orientations,  namely  45  and  22.5  degrees  with  respect  to  the 
east-west  (EW)  principal  direction.  As  shown  in  Fig.  2.11,  the  acceler¬ 
ometers  are  numbered  as  23  and  25  for  the  45-degree  inclination  and  24 
and  26  for  the  22.5-degree  inclination.  A  drawing  for  the  design  of 
the  blocks,  which  are  of  the  same  size  as  the  3-D  blocks,  is  shown  in 
Fig.  2.13.  With  these  two  accelerometer  blocks,  shear  waves  polarized 
at  45  and  22.5  degrees  to  the  east-west  axis  and  propagating  along  the 
vertical  principal  axis  could  be  monitored.  However,  the  study  of  shear 
waves  is  not  covered  in  this  report. 

2.5.4  Charge  Amplifiers 

Since  only  nine  Endevco  charge  amplifiers,  model  2735,  were 
purchased  earlier,  the  29  accelerometers  were  switched  among  them 
before  the  signals  were  displayed  on  the  two  digital  oscilloscopes. 


a)  TOP  VIEW  b)  BOTTOM  VIEW 


Note:  All  dimensions  are  in  inches 


c)  VIEW  A-A 


Fig.  2.13  -  Accelerometer  Block  Used  to  Monitor  Oblique  Shear  Waves 
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Therefore,  these  nine  charge  amplifiers  were  calibrated  to  make  sure 
that  they  provided  the  same  system  gain.  The  charge  amplifiers  were 
calibrated  for  the  full-scale  sensitivity  test  following  the  procedure 
described  In  the  Instruction  manual.  A  schematic  drawing  for  the  elec¬ 
trical  set-up  Is  shown  Fig.  2.14.  The  full-scale  test  was  conducted  for 
a  combination  of  sensitivity  dial  readings  and  full-scale  ranges;  the 
sensitivity  dial  readings  were  1,  3  or  10  pc/g  while  the  full  sca¬ 
le-range  switch  was  set  at  1,  3,  10,  30,  100,  300,  1000  or  3000.  At  all 
times,  the  product  of  the  setting  of  the  sensitivity  dial  and  the 
full-scale  range  switch  was  equal  to  the  number  of  plcocoulombs  (pc) 
which  will  produce  full-scale  output.  The  results  were  satisfactory 
with  minor  adjustment  of  potentiometer  occasionally  necessary.  The 
system  gain  had  an  accuracy  of  better  than  ±1.5  percent  on  all 
full-scale  ranges  which  means  that  the  output  from  any  two  accelerome¬ 
ters  should  be  within  1.5  percent  for  the  same  Input. 

2.5.5  Stress  Cells 

Three  total  stress  cells,  model  TE-9010,  and  a  control  unit, 
model  C-9001,  were  purchased  earlier  from  the  Terra  Technology  at  Red¬ 
mond,  Washington.  (See  section  3.6  of  Kopperman  et  al ,  1982  for  a 
description  of  the  stress  cells.)  Each  stress  cell  was  calibrated 
before  It  was  placed  In  the  sand.  The  set-up  for  calibrating  the 


stress  cells  was  basically  the  same  as  that  used  for  calibrating  the 
strain  gages  on  the  excitation  ports  (Section  2.4).  However,  this  time 
a  larger  capacity  load  cell  (1000-1b  (4450  N))  was  used  to  monitor  the 
applied  load.  In  an  attempt  to  have  a  uniform  pressure  over  the  faces 
of  the  gages,  each  stress  cell  was  covered  with  rubber  on  both  faces 
before  It  was  placed  between  the  top  and  bottom  6-1 n.  (15.2  cm)  square 
platens.  The  platens  had  the  same  size  as  the  plane  dimensions  of  the 
stress  cell.  The  loading  frame  was  then  used  to  load  the  stress  cell  to 
known  pressures  and  readings  from  the  control  unit  of  the  stress  cells 
were  recorded.  Calibration  curves  obtained  are  shown  In  Figs.  2.15  to 
2.17.  The  curves  are  essentially  straight  lines  above  a  pressure  of  5 
psl  (34  kPa)  with  values  of  correlation  coefficient  varying  from  0.99 
to  1.00.  There  Is  no  hysteresis  effect  upon  load  cycling. 

2.5.6  Strain  Sensors 

The  strain  sensors  were  calibrated  using  the  calibration 
fixture  and  control  unit  purchased  from  the  manufacturer.  The  strain 
sensors  are  the  same  as  those  used  before,  namely  model  4000  series. 
Bison  soil  strain  gages  (Section  3.7,  Kopperman  et  al,  1982).  There  are 
four  pairs  of  2-1n.  (5.1  cm)  diameter  strain  sensors  and  four  pairs  of 
4-1n.  (10.2  cm)  diameter  strain  sensors.  Dial  calibration  curves  were 
generated  following  procedures  recommended  In  the  manufacturer's  manu- 
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al  for  each  pair  of  strain  sensors.  With  an  Initial  null  amplitude 
reading  from  the  control  unit  (prior  to  any  strain),  the  calibration 
factor  can  be  obtained  from  the  calibration  curve.  Then,  strains  can  be 
obtained  by  multiplying  this  factor  by  the  difference  In  amplitude 
readings  due  to  soil  movement.  One  of  the  calibration  curves  obtained 
is  shown  In  Fig.  2.18,  and  the  remainder  are  Included  in  Appendix  A. 

As  an  example  of  use  of  the  calibration  curves,  assume  the 
Initial  null  amplitude  for  the  strain  sensor  SN-1  Is  600,  the  corre¬ 
sponding  calibration  factor  Is  0.0505  percent  per  unit  amplitude 
change  as  shown  In  Fig.  2.18.  Now,  assume  the  null  amplitude  for  the 
strain  sensor  SN-1  Is  598  after  a  pressure  of  10  psl  (68.9  kPa)  Is 
applied.  Therefore,  the  change  In  amplitude  reading  Is  -2,  which  multi¬ 
plied  by  the  calibration  factor  (0.0505)  gives  the  resulting  strain  as 
-O.iOl  percent  (negative  for  compression). 

2.6  SUMMARY 

The  same  7-ft  (2.1  m)  cubical  loading  device  used  by  Kopper- 
man  et  al  (1982)  and  Knox  et  al  (1982)  was  used  In  the  present  series  of 
tests,  except  that  the  following  additions  and  Improvements  were  made 
to  the  testing  facilities:  1.  a  bln  for  storing  dry  sand  was  con¬ 
structed,  2.  a  new  ralner  for  placing  the  sand  was  built,  3.  strain 
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gages  were  mounted  on  each  excitation  port  to  control  the  pressure 
applied  at  the  ports,  4.  accelerometers,  stress  cells  and  strain  cells 
were  calibrated,  and  5.  additional  accelerometers  were  Included  to 
provide  more  complete  data  for  an  attenuation  study  and  to  monitor 
shear  waves  polarized  obliquely.  After  the  above  items  were  completed, 
the  sand  specimen  was  built  and  a  series  of  tests  was  commenced  as  dis¬ 
cussed  in  the  following  chapters. 


CHAPTER  THREE 


TESTING  PROGRAM  AND  PROCEDURES 

3.1  INTRODUCTION 

The  main  purposes  of  this  project  were  to  Improve  the  test 
facility  as  discussed  In  Chapter  Two  and  then  to  perform  a  series  of 
tests  on  the  same  sand  used  by  Knox  et  al  (1982)  and  Kopperman  et  al 
(1982). This  series  of  tests  would  Include  all  of  tests  performed  earli¬ 
er  as  well  as  additional  tests  under  anisotropic  states  of  stress. 
These  tests  would  then  be  used  to  check  the  validity  of  the  conclusions 
reported  by  Kopperman  et  al  (1982)  as  well  as  to  extend  the  earlier 
work. 

The  first  step  In  the  test  program  was  to  empty  the  sand  from 
the  trlaxlal  device.  (However,  before  the  device  was  emptied,  Knox  and 
Kopperman  performed  some  tests  with  the  writer  and  S.  Lee  before  they 
left  so  that  the  basic  technique  In  performing  the  tests  could  be  dem¬ 
onstrated.)  Once  the  device  was  emptied,  Improvements  and  calibration 
of  the  Instruments  as  discussed  In  Chapter  Two  were  pursued.  Upon  com¬ 
pletion  of  this  work,  the  device  was  filled  following  the  procedures 
reported  by  Kopperman  et  al  (1982)  with  the  same  sand  that  was  orglnal- 
ly  In  the  device.  Dynamic  testing  was  then  performed  under  the  states 


of  Isotropic  confinement  (5i  =  02  =  53),  biaxial  confinement  (oj  >  o2  =  o3 
or  oj  =  o2  >  53)  ,  and  trlaxlal  confinement  [ai  >  a2  >  o3).  After  each 
sequence  of  testing,  data  were  reduced  and  analyzed  before  performing 
the  next  sequence  of  testing  so  that  additional  Insight  and  better 
understanding  and  planning  for  the  next  sequence  could  be  obtained. 

Kopperman  et  ai  (1982)  and  Hoar  and  Stokoe  (1978)  have  shown 
that  measurements  of  Interval  velocities  (velocities  determined  from 
travel  times  between  receivers)  represent  more  correctly  the  material 
than  direct  velocities  (velocities  determined  from  travel  times 
between  the  source  and  receivers).  Therefore,  only  interval  velocity 
measurements  were  made,  and  the  triggering  system  used  by  Knox  and  Kop¬ 
perman  was  unnecessary  and,  hence,  deleted. 

3.2  EMPTYING  OF  THE  TRIAXIAL  OEVICE 

The  emptying  process  commenced  near  the  start  of  February, 
1982  and  took  about  three  months  to  finish.  The  top  cover  of  the  trlax¬ 
lal  device  was  first  removed  with  the  help  of  an  overhead  crane.  Water 
was  then  bled  from  the  top  membrane  after  which  the  membrane  was 
removed.  The  next  step  In  the  process  was  removal  of  sand.  Sand  was 
removed  with  the  aid  of  a  concrete  bucket  which  was  raised  to  the  top 
of  the  trlaxlal  device  using  a  fork-lift  and  sand  was  shovelled  Into 


51 


the  bucket.  It  was  soon  found  that  this  emptying  process  was  very 
labor  intensive,  especially  after  the  upper  two  ft  (61  cm)  of  sand  was 
removed  and  the  remainder  had  to  be  emptied  with  a  shovelling  action 
that  required  the  shovel  to  be  raised  over  the  shoulder  in  order  to 
dump  sand  into  the  bucket.  The  first  half  of  the  sand  in  the  tri axial 
device  was  stored  in  the  new  storage  bln  (Section  2.2),  and  the  remain¬ 
der  was  stored  in  an  open-air  bin  outside  the  Laboratory. 

Once  the  sand  was  lowered  to  the  level  of  the  embedded 
instrumentation,  the  shovelling  action  became  very  slow  and  tedious. 
The  instrumentation  was  carefully  uncovered,  and  the  locations  of  all 
instruments  were  measured  and  found  to  agree  closely  with  the  locations 
reported  by  Kopperman  et  al  (1982).  These  measurements  were  performed 
as  a  check  on  the  earlier  work  and  to  make  sure  no  mass  movement  of  sand 
had  unknowingly  occurred.  After  all  Instruments  were  removed  and  the 
trlaxial  device  was  completely  emptied,  efforts  were  directed  to  cali¬ 
bration  of  the  Instrumentation  (Sections  2.4  and  2.5). 

Emptying  the  trlaxial  device  took  three  months  because  of: 
1.  careful  checking  of  Instrumentation  locations,  and  2.  Injury  of  S. 
Lee  during  the  emptying  process. 


3.3  SAMPLE  CONSTRUCTION 


A  new  test  specimen  was  constructed  following  the  procedures 
reported  by  Kopperman  et  al  (1982).  The  Inside  of  the  trlaxlal  device 
and  membranes  were  first  cleaned.  One  membrane  was  then  hung  on  the 
north  and  south  walls.  Two  plastic  sheets  with  WD-40  oil  in  between 
were  placed  on  the  bottom  and  four  side  walls  inside  the  triaxial 
device.  Filling  of  sand  Into  the  device  then  started  in  early  August, 
1982. 


The  sand  In  the  new  storage  bln  was  used  first  to  fill  the 
trlaxlal  device.  With  the  small  trap  door  near  the  base  of  the  bin, 
sand  was  allowed  to  flow  out  of  the  bln  and  into  a  concrete  bucket 
which  was  placed  next  to  the  trap  door.  It  took  less  than  ten  minutes 
to  fill  the  bucket  which  was  then  moved  inside  the  Ferguson  Laboratory 
with  a  fork-lift.  An  overhead  crane  was  used  to  lift  the  bucket  so 
that  the  sand  could  be  unloaded  Into  the  new  rainer  which  was  located 
on  the  reinforced  collar  above  the  trlaxlal  device  (Section  2.3).  Once 
the  rainer  was  full  (one  bucket  of  sand  was  enough  to  fill  the  rainer), 
the  trap  doors  under  the  rainer  were  opened  by  about  one  in.  (2.5  cm), 
and  the  rainer  was  moved  from  side  to  side  In  the  east-west  direction 
to  rain  the  sand  Into  the  trlaxlal  device.  It  took  about  twenty  min- 
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utes  to  empty  a  ralner  full  of  sand.  The  raining  process  proceeded  sat¬ 
isfactory  with  the  new  ralner. 

After  all  sand  In  the  storage  bln  had  been  used,  sand  from 
the  open-air  bln  was  used.  Dry  sand  was  shovelled  Into  the  bucket,  and 
the  same  filling  process  followed.  However,  since  there  were  several 
heavy  rainfalls  during  this  time,  sand  In  this  storage  bln  became  wet 
and  had  to  be  spread  out  on  the  concrete  floor  and  dried  In  the  sun 
before  It  could  be  used.  This  drying  operation  slowed  the  filling  pro¬ 
cess  significantly  and  also  made  It  laborious.  (It  Is  recommended  that 
a  second  storage  bln  be  constructed  so  that  all  sand  required  to  fill 
the  trlaxlal  device  can  be  stored  In  dry  areas.)  Moisture  measurements 
showed  that  the  sand  placed  In  the  trlaxlal  device  had  a  moisture  con¬ 
tent  less  than  0.05  percent. 

The  filling  process  was  halted  at  the  eight  elevations  where 
Instruments  were  to  be  placed  as  shown  In  Fig.  2.11.  However,  most 
instruments  were  placed  at  the  mid-depth  plane  as  detailed  in  Fig.  3.1. 
Pictures  of  the  Instruments  are  shown  In  Figs.  3.2  and  3.3.  Since  the 
overhead  crane  was  often  not  available  because  of  other  commlttements 


at  the  Laboratory,  a  wooden  platform  was  hung  Inside  the  trlaxlal 
device  with  ropes  tied  to  the  four  castors  of  the  ralner  as  shown  In 
Fig.  3.4.  This  was  done  so  that  one  could  climb  (carefully)  Into  the 
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Legend:  1.  □  2-D  or  3-D  Accelerometers 

2.  SS-1  to  SS-3  stress  cells  (which  are  actually 
6  in.  above  the  mid-depth  plane) 

3.  SN-1  to  SN-4:  2  in.  strain  sensors 

4.  SN-5  to  SN-6:  4  in.  strain  sensors 

Fig.  3.1  -  Instrumentation  at  Mid-Depth  Plane  in  Triaxial 
Device 
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Legend:  1.  IS  to  13W:  accelerometer  levels;  letter  following  accelero¬ 
meter  number  indicates  direction  of  sensivity  (see  Fig.  2.11) 

2.  SN-1  to  SN-4:  2-in.  strain  sensors 

3.  SN-6  :  4-in.  strain  sensors 

Fig.  3.2  -  Instruments  Being  Placed  at  the  Mid-Depth  Plane 


trlaxlal  device  and  use  the  platform  as  a  working  area  without  disturb¬ 
ing  the  sand.  This  platform  was  used  to  place  all  embedded  Instrumenta¬ 
tion. 

The  filling  process  was  also  stopped  at  four  elevations  whe¬ 
re  density  measurements  were  taken.  The  same  density  devices  used  by 
Knox  and  Kopperman  were  used.  These  devices  are  three  cylindrical  con¬ 
tainers  made  of  0.25-1n.  (0.64  cm)  thick  plexiglass.  The  containers 
have  an  Internal  diameter  of  6  In.  (15.2  cm)  and  a  height  of  3  In.  (7.6 
cm),  which  results  In  a  volume  of  8.64  In.1  (0.05  ft*  or  142  cm1).  The 
top  edge  of  the  cylinder  Is  machined  to  an  angle  of  45  degrees  sloping 
out  to  minimize  sand  bouncing  Into  the  container  upon  striking  the  con¬ 
tainer  during  the  raining  process,  which  might  alter  the  compactness  of 
the  density  specimen.  The  elevations  and  locations  for  density  meas¬ 
urements  are  shown  In  Fig.  3.5,  and  the  results  are  listed  In  Table 
3.1.  The  average  density  Is  101.8  pcf  (1629  kg/m1),  and  the  standard 
deviation  Is  2.0  pcf  (32  kg/m1),  less  than  2  percent  of  the  mean  value. 
Compared  to  values  obtained  by  Kopperman  et  al,  1982  (  mean  value  and 
standard  devltalon  of  96.6  pcf  and  1.9  pcf,  respectively),  the  sample 
Is  denser  than  before,  with  density  5.4  percent  higher.  In  terms  of 
average  void  ratio,  the  new  sample  has  a  value  of  0.64  and  Is  12  per¬ 
cent  lower  than  the  value  of  0.73  reported  by  Kopperman  et  al.  The  val¬ 
ues  of  standard  deviation  show  the  same  amount  of  experimental  scatter 


(a)  Plan  View  of  Measurement  Locations 


(b)  Elevations  Where  Raining  of  Sand 
was  Halted  for  Measurements 

Fig.  3.5  -  Locations  of  Density  Measurements  Performed 
During  Filling  of  Triaxial  Device 
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for  both  samples.  The  new  ralner  apparently  builds  a  denser  sample  but 
not  necessarily  a  more  homogeneous  sample  than  before. 


3.4  DYNAMIC  TESTING 

Construction  of  the  sample  was  finished  In  early  October, 
1982.  A  series  of  dynamic  tests  was  then  performed  for  Isotropic  con¬ 
finement,  followed  by  a  series  of  tests  with  biaxial  confinement. 
Unfortunately  the  top  membrane  started  to  leak  In  early  January,  1983. 
This  forced  testing  to  be  halted.  The  top  membrane  was  removed,  and 
locations  of  leakage  were  found.  Due  to  the  special  construction  of  the 
membrane,  the  repair  cost  was  unreasonable  compared  to  the  cost  of  a 
new  membrane.  Therefore,  a  new  membrane  was  ordered  from  Goodyear  Aero¬ 
space  Corporation  at  Rockmart,  Georgia,  the  same  company  that  built  the 
old  membranes.  In  April,  1983,  the  new  membrane  was  dellevered.  Testing 
was  resumed  In  May,  1983,  and  was  completed  In  late  May,  1983. 

The  loading  history  for  all  series  of  tests  Is  summarized  In 
Table  3.2.  The  history  Is  also  shown  In  Fig.  3.6  for  the  Isotropic  and 
biaxial  testing  series  before  the  membrane  leaked  and  Fig.  3.7  for  the 
biaxial  and  trlaxlal  testing  series  after  the  new  membrane  was 


Installed. 


TABLE  3.2 

Loading  Pressure  Sequences 


Horizontal 


Test 

No. 

Date  of  Test 

Vertical  Effective  Stress 
o TB  (psi)* 

Effective 
oNS  (psi)* 

Stresses 
5EW  (psi)* 

1 

10/08/82 

10 

10 

10 

2 

10/13/82 

15 

15 

15 

3 

10/16/82 

20 

20 

20 

4 

10/16/82 

30 

30 

30 

5 

10/20/82 

40 

40 

40 

6 

10/23/82 

30 

30 

30 

7 

10/23/82 

20 

20 

20 

8 

10/23/82 

15 

15 

15 

9 

10/23/82 

10 

10 

10 

10 

11/08/82 

15 

15 

15 

11 

11/10/82 

20 

15 

15 

12 

11/12/82 

30 

15 

15 

13 

11/13/82 

40 

15 

15 

14 

11/13/82 

40 

20 

20 

15 

11/15/82 

40 

30 

30 

16 

11/28/82 

40 

40 

40 

17 

12/17/82 

40 

30 

30 

18 

12/17/82 

40 

20 

20 

19 

12/17/82 

40 

15 

15 

20 

12/18/82 

30 

15 

15 

21 

12/18/82 

20 

15 

15 

22 

12/18/82 

15 

15 

15 

23 

12/18/82 

20 

20 

20 

24 

12/20/82 

30 

20 

20 

25 

12/20/82 

40 

20 

20 

26 

12/20/82 

15 

20 

20 

27 

12/21/82 

10 

20 

20 

28 

12/21/82 

20 

10 

20 

29 

12/21/82 

20 

15 

20 

30 

12/22/82 

20 

30 

20 

31 

12/22/82 

20 

40 

20 

Date  of  Test 


TABLE  3.2 

Loading  Pressure  Sequences 
(Continued) 


Vertical  Effective  Stress 
* 


Horizontal 
Effective  Stresses 
* 


12/22/82 

30 

01/05/83 

40 

01/06/83 

20 

01/06/83 

20 

01/06/83 

20 

01/07/83 

20 

01/07/83 

15 

01/07/83 

15 

01/08/83 

15 

01/08/83 

15 

01/08/83 

20 

01/10/83 

30 

01/10/83 

20 

01/10/83 

15 

01/11/83 

15 

01/11/83 

15 

01/11/83 

15 

01/11/83 

15 

01/12/83 

15 

01/12/83 

20 

01/12/83 

30 

05/19/83 

30 

05/20/83 

40 

05/20/83 

30 

05/20/83 

20 

05/21/83 

15 

05/23/83 

15 

05/23/83 

15 

05/24/83 

15 

<v.*, 
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Table  3.2 


Loading  Pressure  Sequences 
(Continued) 


Horizontal 

Vertical  Effective  Stress  Effective  Stresses 
Date  of  Test  5TB  (psi)*  oNS  (psi)*  5EW  (psi)* 


05/24/83 

40 

25 

20 

05/24/83 

40 

15 

25 

05/24/83 

40 

15 

30 

05/24/83 

40 

15 

35 

05/25/83 

28 

28 

28 

05/25/83 

32 

28 

24 

05/25/83 

36 

28 

20 

05/25/83 

40 

28 

16 

05/26/83 

28 

28 

28 

05/26/83 

32 

31 

24 

05/26/83 

36 

34 

20 

05/26/83 

40 

37 

16 

TB  *  Vertical  (top-bottom)  effective  stress 
'NS  *  Horizontal  (north-south)  effective  stress 
EW  *  Horizontal  (east-west)  effective  stress 


Fig.  3.7  -  Loading  History  in  Each  Principal  Stress 
Direction  After  New  Membrane  Installation 
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3.5  COMPRESSION  WAVE  ANALYSIS 

Waveforms  of  the  compression  wave  were  monitored  under  each 
confining  pressure  and  were  recorded  on  magnetic  disks  using  two  dig¬ 
ital  oscilloscopes.  These  records  were  brought  back  to  the  main  campus 
for  analyses  from  which  propagation  velocities,  frequencies,  and 
strain  amplitudes  were  determined  with  procedures  detailed  In  the  fol¬ 
lowing  sections. 

3.5.1  Determination  of  Propagation  Velocity 

With  known  distances  between  accelerometers  (which  were  mea¬ 
sured  from  the  center  of  one  accelerometer  to  the  center  of  next  accel¬ 
erometer  during  sample  construction)  and  measured  travel  times  of  the 
waves,  wave  propagation  velocities  were  simply  calculated  by  dividing 
the  distances  by  the  corresponding  travel  times.  The  distances  used  In 
the  analysis  are  summarized  In  Table  3.3.  The  waveforms  recorded  on  the 
magnetic  disks  were  recalled  using  the  same  digital  oscilloscopes  used 
In  collecting  the  waves,  and  travel  times  were  estimated  from  first 
arrivals  of  the  waves  at  each  accelerometer.  A  plotting  of  typical 
waveforms  Illustrating  this  procedure  Is  shown  In  Fig.  3.8. 
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TABLE  3.3 

Distances  Between  Accelerometers  Inside  the  Triaxlal  Device 


Accelerometer  Labels* 
From  To 


IS 

4S 

1.13 

4S 

6S 

0.99 

6S 

9S 

0.99 

9S 

12S 

2.00 

16W 

10W 

2.07 

10W 

19W 

2.06 

20V 

8V 

2.00 

8V 

27V 

2.00 

2W 

5W 

1.13 

5W 

7W 

0.99 

7W 

10W 

0.99 

10W 

13W 

2.00 

3V 

8V 

1.98 

8V 

11V 

2.00 

15V 

9S 

2.07 

9S 

18S 

2.06 

14V 

8V 

2.07 

8V 

17V 

2.06 

21S 

9S 

2.00 

9S 

28S 

2.00 

22W 

10W 

2.00 

10W 

29W 

2.00 

23 

25 

2.00 

24 

26 

2.00 

*See  Fig.  2.11  for  Accelerometer  Locations 
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3.5.2  Determination  of  Predominant  Frequency 

The  predominant  frequency  of  the  compression  wave  was  deter¬ 
mined  from  two  fractional  values  of  the  period,  0.25T  and  0.50T,  meas¬ 
ured  from  each  accelerometer  record  as  shown  In  Fig.  3.9.  These 
frequencies  were  then  used  to  estimate  strain  amplitudes  and  the  wave¬ 
lengths  as  discussed  In  Sections  3.5.4  and  3.5.5,  respectively.  Only 
rough  approximations  of  the  values  were  of  concern.  Hence,  only  four 
sets  of  records  were  used  In  the  analysis  :  two  sets  for  the  isotropic 
state  of  confinement,  one  set  for  the  biaxial  state  of  confinement  and 
one  set  for  the  triaxlal  state  of  confinement.  Each  set  of  records 
contained  all  the  compression  waves  propagating  along  each  principal 
axis  under  the  particular  stress  state. 

Values  of  predominant  frequency  determined  from  0.25T  varied 
from  500  Hz  to  2632  Hz  with  an  average  of  1287  Hz  and  standard  devi¬ 
ation  of  426  Hz.  Values  of  predominant  frequency  determined  from  0.50T 
varied  from  1205  Hz  to  5882  Hz  with  an  average  of  2400  Hz  and  standard 
deviation  of  1044  Hz.  Frequencies  determined  from  0.50T  were  much  high¬ 
er  than  those  from  0.25T  due  to  the  phenomenon  of  a  slow  rise  time  of 
the  initial  pulse.  The  average  frequency  for  all  measured  values  from 
both  0.25T  and  0.50T  was  1830  Hz. 


3.5.3  Determination  of  Peak  Acceleration 

Two  amplitudes  of  each  wave  form  were  determined  from  the 
same  four  sets  of  records.  The  amplitudes  determined  were  the 
first-peak  (or  trough)  amplitude,  Al,  and  the  first  peak-to-peak 
amplitude,  A 2,  as  shown  In  Fig.  3.9. 

The  charge  amplifiers  used  with  the  accelerometers  were  cal¬ 
ibrated  so  that  the  product  of  the  sensitivity  dial  and  the  full-scale 
switch  was  equal  to  the  number  of  plcocoulombs  (pc)  which  will  produce 
the  full-scale  output  of  20V  peak-to-peak.  As  an  example,  let  the  pro¬ 
duct  be  100  pc  and  the  sensitivity  of  a  typical  accelerometer  be  100 
pc/g.  Therefore,  an  acceleration  of  one  g  (32.2  ft/s*  or  9.8  m/s*) 
peak-to-peak  would  produce  an  output  of  20V.  With  measured  values  of 
Al  and  A 2  In  volts,  peak  and  peak-to-peak  accelerations  of  each  wave 
were  determined.  Values  of  peak  acceleration  ranged  from  0.007g  to 
0.564g  (0.23  to  18.0  ft/s*  or  0.069  to  5.5  m/s*)  with  an  average  of 
0.124g  (3.99  ft/s*  or  1.22  m/s*).  The  average  of  peak-to-peak  acceler¬ 
ations  was  0.275g  (8.86  ft/s*  or  2.70  m/s*)  and  was  about  twice  the 
peak  acceleration. 
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3.5.4  Determination  of  Strain  Amplitude 


It  Is  well  known  that  the  linear  strain  amplitude,  e,  for  a 
plane  wave  can  be  determined  from  the  particle  velocity, Qp , and  propa¬ 


gation  velocity,!/  , as  follows  : 

e  =  ^ 

VP 


(3.1) 


up  =  (2irf)A 

where:  f  =  wave  frequency,  and 

A  =  peak  particle  amplitude. 


(3.2) 


Using  a  harmonic  approximation,  the  peak  acceleration,  z, 
can  be  related  to  peak  amplitude  as  follows: 


Z  =  (2irf)2A 


(3.3) 


which  means: 


0  =  -JL. 

P  2irf 


(3.4) 


From  Eqs.  3.1  and  3.4,  the  relationship  between  strain  amplitude, 
acceleration,  frequency,  and  compression  wave  velocity  for  a  plane 
wave  Is: 


*  =  (2wf)Vp 


(3.5) 
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The  sources  In  the  trlaxial  device  tend  to  generate  spheri¬ 
cal  waves  while  Eq.  3.5  is  derived  for  plane  waves  with  harmonic  exci¬ 
tation.  Therefore,  values  of  strain  amplitude  computed  by  the  equation 
can  only  be  used  as  an  approximation  to  reflect  the  order  of  magnitude 
of  the  strains  in  the  sand  in  the  trlaxial  device.  To  estimate  the 
maximum  strain  amplitude,  the  smallest  values  of  frequency  and  com¬ 
pression  wave  velocity  together  with  the  largest  peak  acceleration 
amplitude  gives  a  strain  amplitude  of  about  0.0005  percent.  Since  this 
is  less  than  0.001  percent,  testing  may  be  considered  to  be  low  ampli¬ 
tude  and  the  effect  of  strain  amplitude  can  be  ignored. 


3.5.5  Determination  of  Wavelengths 


With  known  values  of  compression  wave  velocity,  Vp  ,  and  fre¬ 
quency,  f,  the  wavelength,  X,  can  be  computed  from  the  following 
equation: 


X= 


(3.6) 


Using  the  average  minimum  velocity  of  1340  ft/s  (409  m/s)  and  the  aver¬ 
age  maximum  velocity  of  1980  ft/s  (604  m/s)  together  with  an  average 


frequency  of  1830  Hz,  the  range  In  wavelengths  was  0.73  ft  to  1.08  ft 
(0.22  m  to  0.33  m)  In  the  sand  In  the  trlaxlal  device. 

As  discussed  In  Section  6.2.5  by  Kopperman  et  al  (1982),  the 
wavelengths  must  be  short  enough  so  that  the  full  wave  passes  by  any 
receiver  before  any  reflected  wave  reaches  the  receiver  In  order  to 
obtain  waveforms  that  are  easily  Interpretable.  This  critical  dis¬ 
tance  Is  twice  the  distance  between  the  receiver  furthest  from  the 
source  and  the  wall  opposite  the  source.  In  this  series  of  tests,  the 
accelerometers  were  positioned  such  that  this  distance  Is  about  3  ft 
(0.92  m)  which  Is  more  than  twice  the  longest  wavelength.  Therefore, 
the  wavelengths  are  short  enough  to  be  read  without  Interference  from 
reflections. 

Conversely,  the  wavelength  should  be  long  enough  so  that  the 
accelerometer  block  properly  tracks  particle  motions  In  the  sand.  If 
the  wavelength  Is  too  small,  the  block  acts  as  a  boundary  rather  than 
as  a  part  of  the  sand  specimen.  The  accelerometer  blocks  embedded  In 
the  sand  are  1.75  In.  (44  mm)  square.  Since  the  minimum  wavelength  Is 
0.73  ft  (8.8  In.  or  220  mm)  which  Is  five  times  the  dimension  of  the 
block,  wavelengths  may  be  considered  as  long  enough  to  have  waves  that 
are  easy  to  Interpret  (Suddhlprakarn  et  al ,  1983). 


3.6  SUMMARY 


The  trlaxlal  device  used  by  Knox  et  al  (1982)  and  Kopperman 
et  al  (1982)  was  first  emptied.  Improvements  and  calibration  of  the 
Instruments  as  discussed  In  Chapter  Two  were  pursued.  Then,  a  sample  of 
the  same  sand  was  constructed  using  the  new  raining  device  (Section 
2.3).  The  filling  process  was  halted  at  elevations  where  Instruments 
were  to  be  placed  and  density  measurements  were  taken.  The  locations  of 
Instruments  are  shown  In  Figs.  2.11  and  3.1.  The  average  density  of  the 
sample  Is  101.8  pcf  (1629  kg/m1). 

Once  the  sample  was  constructed,  series  of  dynamic  tests 
were  performed  for  Isotropic  confinement,  biaxial  confinement  and  trl¬ 
axlal  confinement.  The  loading  history  for  all  series  of  tests  Is  sum¬ 
marized  In  Table  3.2  and  shown  In  Figs.  3.6  and  3.7.  Analysis  for 
waveforms  of  the  compression  wave  collected  under  each  test  Is  pre¬ 
sented  In  Section  3.5.  Propagation  velocity  of  each  P-wave  Is  deter¬ 
mined  and  relationship  between  P-wave  velocity  and  states  of 
confinement  Is  studied  and  presented  In  Chapters  4,  6  and  7. 


The  first  step  in  the  study  of  the  effect  of  confining  pres¬ 
sure  on  compression  wave  velocity  was  to  perform  tests  with  isotropic 
confinement.  This  state  of  confinement  is  the  simplest  state  that  can 
be  applied  with  the  triaxlal  device,  and  the  effect  of  stress  state  on 
velocity  is  the  easiest  to  analyze.  Moreover,  this  state  of  confinement 
eliminates  stress  as  a  contributing  factor  for  different  values  of 
P-wave  velocity  measured  along  each  of  the  principal  axes  and  thus  per¬ 
mits  evaluation  of  structural  anisotropy  (also  called  Inherent  aniso¬ 
tropy)  . 

Compression  wave  velocities  were  measured  at  various  iso¬ 
tropic  states  of  stress  ranging  from  10  psi  (69  kPa)  to  40  p$1  (276 
kPa).  To  obtain  Isotropic  loading,  confining  pressures  were  applied  so 
that  stresses  along  all  three  principal  axes  of  the  triaxlal  device 
were  the  same  (a1=a2=53).  However,  In  reality,  true  Isotropic  con¬ 
finement  could  not  be  achieved  at  all  points  within  the  triaxlal  device 
because  of  small  Inherent  pressure  gradients.  In  the  vertical  direc¬ 
tion,  there  was  a  slight  Increase  in  stress  with  depth  caused  by  the 
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weight  of  sand.  The  change  In  vertical  stress  from  top  to  bottom  of  the 
device  was  less  than  5  psl  (34.5  kPa),  or  less  than  ±2.5  psi  (17  kPa) 
from  the  value  applied  at  the  center  of  the  device.  For  stresses 
applied  In  the  horizontal  directions,  the  hydraulic  head  of  water  In 
the  membranes  caused  a  change  in  stress  from  top  to  bottom  of  the 
device  of  about  3  psi  (20.7  kPa),  or  less  than  ±1.5  psi  (10.3  kPa)  from 
the  value  applied  at  the  mld-helght  of  the  device.  The  largest  verti¬ 
cal  distance  away  from  the  mid-height  of  the  specimen  to  any  accelerom¬ 
eter  was  2  ft  (0.61  m)  as  shown  In  Fig.  2.11.  Therefore,  variations  in 
pressures  in  the  area  over  which  data  were  collected  were  ±1.4  psi  (9.6 
kPa)  and  ±0.86  psi  (5.9  kPa)  in  the  vertical  and  horizontal  directions, 
respectively.  Within  the  range  of  Isotropic  pressures,  namely  10  to  40 
psi  (69  to  276  kPa),  variations  In  pressures  over  the  area  of  interest 
were  3.5  to  14  percent  and  2.2  to  8.6  percent  In  the  vertical  and  hori¬ 
zontal  directions,  respectively.  These  pressure  gradients  were  small 
and  were  considered  negligible. 

Values  of  all  pressures  reported  in  this  study  represent  the 
values  of  pressures  applied  at  the  center  of  the  trlaxlal  device.  To 
account  for  the  effects  of  overburden  pressure,  hydraulic  gradient, 
and  different  piezometric  elevations  between  the  membranes  and  pres¬ 
sure  panel  board,  the  following  corrections  were  applied  to  the  pres¬ 
sure  gage  readings  : 


vertical  correction  :  pg  =  pc  -  0.21  psl 

horizontal  correction  :  p  =  P  +  0.72  psl 

G  C 

where  :  Pg  =  reading  on  pressure  gage  at  panel  board,  in  psi,  and 
=  pressure  at  center  of  triaxlal  device,  in  psi. 

These  corrections  were  derived  for  a  constant  difference  in  elevation 
between  the  top  or  center  of  the  device  and  the  air/water  Interface  at 
the  panel  board,  namely  for  a  10-in.  (25  cm)  high  column  of  water  in 
the  accumulator  on  the  panel  board.  It  was  difficult  to  keep  the  water 
column  In  the  accumulator  at  a  constant  elevation  throughout  measure¬ 
ments  at  a  particular  set  of  confining  pressures.  A  change  in  elevation 
of  the  water  level  in  the  accumulator  of  one  in.  (2.5  cm)  would  result 
in  about  a  0.04  psi  (0.25  kPa)  change  In  pressure,  and  a  change  of  the 
entire  18-in.  (45.7  cm)  high  accumulator  would  result  in  only  a  0.65 
psi  (4.48  kPa)  change  In  pressure.  Therefore,  with  the  effort  that  was 
made  to  maintain  the  water  level  at  an  approximately  constant  ele¬ 
vation,  any  set  of  pressures  can  be  considered  constant  (within  0.2  psi 
(1.4  kPa)) . 


The  method  of  data  analysis  is  discussed  In  Section  3.5.  The 
effect  of  stress  history  on  P-wave  velocity  was  studied  first.  Then, 
the  effect  of  structural  (Inherent)  anisotropy  and  variation  of  P-wave 
velocity  with  isotropic  loading  were  studied. 
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4.2  EFFECT  OF  STRESS  HISTORY 

An  Initial  series  of  tests  was  performed  under  isotropic 
loading  with  confining  pressures  Increased  from  10  psi  (69  kPa)  to  40 
psl  (276  kPa).  The  sand  specimen  was  then  unloaded  Isotropically  to  10 
psi  (69  kPa).  At  each  state  of  isotropic  confinement,  compression  wave 
velocities  were  determined  as  the  average  of  two  interval  velocities 
(Section  3.1)  along  each  of  the  principal  axes,  except  for  the 
north-south  axis  along  which  three  interval  velocities  (provided  by 
additional  accelerometers  used  for  an  attenuation  study)  were  used  to 
obtain  the  average  velocities. 

It  is  worthwhile  to  examine  first  the  scatter  in  the  interval 
velocities  along  each  axis.  Along  the  north-south  axis,  the  Intermedi¬ 
ate-interval  velocity  was  always  less  than  the  near-interval  velocity, 
with  the  average  value  of  the  ratio  of  these  two  velocities  being  0.98. 
On  the  other  hand,  the  far-interval  velocity  was  always  greater  than 
the  intermediate-interval  velocity.  The  ratio  of  these  two  velocities 
averaged  1.14.  Similarly,  the  ratio  of  the  far-interval  velocity  to 
the  near-interval  velocity  along  the  vertical  axis  had  an  average  value 
of  1.08.  Along  the  east-west  axis,  the  ratio  of  the  far-interval 
velocity  to  the  near-interval  velocity  averaged  1.005.  All  of  these 
variations  show  the  Inherent  scatter  in  the  system.  Since  the  differ- 
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ences  were  not  large,  average  values  of  the  Interval  velocities  were 
used  In  all  analyses. 


Compression  wave  velocities  determined  from  averages  of 
interval  velocities  were  plotted  against  the  corresponding  confining 
pressure  on  log-log  paper  for  each  of  the  principal  axes  as  shown  In 
Fig.  4.1.  Straight  lines  were  fitted  to  the  data  for  loading  and 
unloading  separately  using  a  least-squares  method.  As  shown  In  Fig. 
4.1,  straight  lines  for  the  unloading  data  have  slightly  flatter  slopes 
than  those  for  the  loading  data  which  results  In  the  values  of  m  being 
greater  upon  loading.  This  effect  Is  most  Important  at  the  lowest 
pressure.  However,  the  largest  variation  In  VpUpon  loading  and  unload¬ 
ing  is  only  6  percent.  The  hysteresis  effect  upon  was,  therefore, 
not  significant,  and  averages  of  values  obtained  from  both  loading  and 
unloading  were  used  In  all  analyses. 

Testing  at  an  Isotropic  state  of  stress  occurred  at  various 
times  throughout  the  complete  series  of  tests  as  summarized  In  Table 
3.2  or  as  shown  In  Figs.  3.6  and  3.7.  In  an  attempt  to  study  further 
the  effect  of  stress  history  on  P-wave  velocity,  velocities  from 
repeated  tests  are  summarized  In  Table  4.1.  Averages  and  standard  devi¬ 
ations  of  velocities  at  each  Isotropic  state  of  stress  were  determined 
for  each  of  the  three  principal  axes.  As  can  be  seen  In  the  table,  com- 
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c)  Propagation  in  Vertical  (TB)  Direction 


Confining  Pressure,  oQ  (psi) 

Fig.  4.1  -  Effect  of  Stress  History  on  Variation  of 
P-wave  Velocity  with  Isotropic  Confining 
Pressure 


'  V  v  < 


TABLE  4.1 


Compression  Wave  Velocities  Measured  at  Different  Times 
Under  Similar  Isotropic  States  of  Stress 


P-wave  Velocities  (fps) 


Confining 

Test* 

HORIZONTAL 

DIRECTIONS 

VERTICAL  DIRECTION 

Pressure 

No. 

(north-south) 

(east-west) 

(top-bottom) 

10  psi 

1 

1434 

1416 

1270 

9 

1503 

1498 

1301 

M**  -  1469 

M  =  1457 

M  =  1286 

SD**  -  49 

SD  -  58 

SD  »  22 

15  psi 

2 

1588 

SifSB 

1413 

8 

1585 

1413 

10 

1624 

MgsMBSm 

1406 

22 

1604 

B mm 

1409 

38 

1606 

1563 

1389 

60 

1607 

1666 

1341 

M  -  1602 

M  -  1569 

M  -  1395 

SD  *  14 

SD  «  49 

SD  -  28 

20  psi 

3 

1695 

1648 

1481 

7 

1758 

1737 

1491 

23 

1697 

1667 

1488 

35 

1798 

1669 

1458 

M  =  1737 

M  -  1680 

M  -  1480 

SD  =  50 

SD  -  39 

SD  -  15 

30  psi 

4 

1875 

1801 

IBS  '  ' 

6 

1881 

1880 

32 

1897 

1891 

53 

1797 

1894 

M  =  1863 

M  -  1867 

SD  -  44 

SD  -  44 

mmm 

40  psi 


5 

16 

33 

34 


1985 

1932 

1688 

1985 

1932 

1688 

2012 

1979 

1665 

1963 

1998 

1669 

M  =  1986 

M  =  I960 

M  -  1678 

SD  -  20 

SD  ■  34 

SD  - _ 12 

pression  wave  velocities  at  each  confining  pressure  are  close.  The 
maximum  percentage  differences  between  any  two  waves  velocities 
propagating  along  the  same  axis  under  the  same  confining  pressure 
are  small,  less  than  8  percent  for  all  cases.  It  should  also  be 
noted  that  the  standard  deviations  are  always  less  than  4.0  percent 
of  the  average  velocities.  Therefore,  the  scatter  of  average 
velocities  can  be  considered  small,  and  the  effect  of  stress 
history  on  P-wave  velocity  Is  negligible. 

4.3  EFFECT  OF  STRUCTURAL  ANISOTROPY 

Structural  anisotropy  Is  the  Inherent  anisotropy  In  the 
soil  skeleton  which  causes  a  difference  In  soil  properties  (Including 
P-wave  velocities)  In  different  directions  under  isotropic  loading 
conditions.  (Soil  properties  will  vary  with  direction  if  the 
stresses  vary,  and  this  Is  referred  to  as  stress-induced  anisotropy.) 
A  certain  amount  of  structural  anisotropy  was  believed  to  have  been 
Induced  In  the  sample  as  a  result  of  the  raining  operation  used  to 
build  the  sample.  Therefore,  it  was  necessary  to  evaluate  the 
magnitude  and  importance  of  the  effect  of  structural  anisotropy  on 
P-wave  velocities  before  the  effect  of  isotropic  confinement  could 
be  studied. 


The  easiest  way  to  study  the  effect  of  structural  anisotropy 
on  P-wave  velocity  Is  to  perform  tests  under  Isotropic  confinement. 
Therefore,  average  P-wave  velocities  summarized  In  Table  4.1  were 

plotted  against  confining  pressure  using  logarithmic  paper.  The 

/ 

resulting  curves  are  shown  In  Fig.  4.2  and  are  represented  by  the  three 
best-fit  straight  lines  that  are  almost  parallel.  The  slopes  of  the 
straight  lines  are  0.22,  0.23 and  0.17  for  velocities  propagating  along 
the  north-south  (NS),  east-west  (EW)  and  vertical  (TB)  directions, 
respectively.  Additionally,  the  curves  fall  Into  two  groups:  1.  the  two 
curves  for  waves  propagating  along  the  horizontal  directions  (NS  and 
EW)  In  which  case  the  velocities  are  within  3.5  percent,  and  2.  the 
curve  for  waves  propagating  along  the  vertical  direction  which  falls 
below  the  other  two  curves.  Therefore,  P-wave  velocities  along  the  NS 
and  EW  directions  can  be  considered  the  same.  As  such,  no  structural 
anisotropy  Is  present  In  the  horizontal  plane  which  can  be  considered  a 
plane  of  Isotropy.  However,  the  ratio  of  velocities  of  P-waves  propa¬ 
gating  along  the  horizontal  directions  (VDU)  to  those  along  the  vertl- 
cal  direction  (Vpy)  varied  from  1.16  at  10  psl  (69  kPa)  to  1.18  at  40 
psl  (276  kPa).  Therefore,  structural  anisotropy  exists  between  the 
horizontal  and  vertical  directions.  This  condition  Is  referred  to  as 
cross-anisotropy.  (See  Section  5.5  for  numerical  evaluation  of  this 
constitutive  model.)  The  average  of  the  ratios  of  Vou  to  VDU  Is 


1.17.  Hence,  the  wave  velocities  may  be  approximated  by  the  following 
relationship  which  Is  Independent  of  stress  (provided  the  stress  level 
Is  within  10  to  40  psl  (69  to  276  kPa)): 

VPH  =  1-17  (Vpy)  (4.1) 

Since  a  pressure  gradient  exists  In  the  trlaxlal  device 
(Section  4.1),  velocities  will  Increase  slightly  with  depth  causing 
some  wave  refraction  (curved  wave  path  for  waves  propagating  along  the 
horizontal  directions).  However,  the  effect  of  a  curved  wave  path  would 
be  small  since  the  receivers  are  closely  spaced  when  compared  to  the 
relative  depth  of  embedment  based  on  confining  pressure  (Haupt,  1973). 
The  differences  In  P-wave  velocities,  therefore,  are  attributable  to 
structural  anisotropy  since  the  same  Instrumentation  and  testing  tech¬ 
nique  were  used  to  determine  P-wave  velocities  propagating  along  each 
principal  axis. 

4.4  EFFECT  OF  ISOTROPIC  CONFINEMENT 

For  the  dry  sand  sample,  confining  pressure  Is  the  same  as 
effective  stress  (pore  water  pressure  Is  zero).  Also,  for  Isotropic 
confinement,  the  confining  pressure  along  each  of  the  three  principal 
axes  Is  the  same  as  the  mean  effective  stress,  5q.  Therefore,  only  one 
velocity  versus  pressure  plot  Is  needed  In  the  analysis.  This  plot  Is 


shown  In  Fig.  4.2  in  which  the  variation  in  P-wave  velocities  with  <jQ 
can  be  represented  by  a  linear  relationship  on  a  log-log  plot  for  all 
three  axes.  The  relation  shown  in  Fig.  4.2  can  be  expressed  in  power 
form  as  used  by  Hardin  and  Richart  (1963): 

Vp  =  (4.2) 

where:  Vp  =  P-wave  velocity,  in  fps, 

aQ  -  isotropic  confining  pressure,  in  psf, 

=  constant  (C^  =  Vp  when  oQ  *  1  psf),  and 
m  «  slope  of  the  straight  line. 

The  slopes  of  the  straight  lines  for  the  NS  and  EW  directions 
have  values  0.22  and  0.23  while  that  for  the  TB  direction  has  a  value 
of  0.17.  These  values  are  slightly  less  than  those  reported  by  Hardin 
and  Richart  (1963).  They  concluded  that  P-wave  velocity  in  sands  varied 
as  in  Eq.  4.2  but  with  a  slope  of  0.25  -  10  percent  within  the  range  of 
2000  psf  to  8000  psf  (96  to  383  kPa).  The  values  of  the  constant  are 
very  close  to  each  other,  namely  308,  280  and  371  for  NS,  EW  and  TB 
directions,  respectively.  (If  the  slope  of  the  straight  line  for  the  TB 
direction  were  also  0.22,  the  value  of  the  constant  back-calculated 
from  P-wave  velocity  at  20  psi  (138  kPa)  would  be  241  which  looks  more 
reasonable  when  compared  to  308  and  280  for  the  horizontal  directions 
for  this  cross-anisotropic  specimen.)  If  P-wave  velocity  and  isotro- 


ic  confining  pressure  are  in  units  of  m/s  and  kPa  respectively,  the 
values  of  the  constant  would  be  183,  172,  and  190  for  NS,  EM  and 
TB  directions,  respectively.  (The  value  of  is  143  in  the  TB  direction 
if  a  slope  of  0.22  is  assumed.) 

With  conditions  of  isotropic  confinement,  differences  between 
velocities  in  the  horizontal  and  vertical  directions  are  due  to 
structural  anisotropy  which  is  seen  to  be  independent  of  confining 
pressure.  Writing  Eq.  4.2  for  velocities  in  the  horizontal  and  vertical 
directions  separately  and  eliminating  the  confining  pressure  from  the 
two  equations,  the  relationship  between  horizontal  and  vertical  veloc¬ 
ities  can  be  expressed  as: 

VpH  «  0.117 ( Vyp)l .323  (4.3) 

It  should  be  noted  that  the  above  equation  is  only  valid  for  the. sand 
sample  used  in  the  present  series  of  tests  since  structural  anisotropy 
can  vary  significantly  with  sample  preparation  (Mul ill’s  et  al ,  1977). 

4.5  SUMMARY 

From  this  study  of  P-wave  velocity  under  isotropic  confine¬ 
ment,  the  following  conclusions  can  be  drawn. 


1.  The  effect  of  stress  history  on  P-wave  velocity  Is  small  and 
can  be  neglected  for  this  sand. 

2.  Structural  anisotropy  affected  the  measured  wave  velocities, 
with  velocities  In  the  two  horizontal  directions  about  17  per¬ 
cent  higher  than  the  velocity  In  the  vertical  direction. 

3.  The  relation  between  P-wave  velocity  and  Isotropic  confining 
pressure  agrees  with  that  proposed  by  Hardin  and  R1 chart  (1963) 
and  can  be  expressed  by  Eq.  4.2.  The  general  form  of  the 
equation  was  unaffected  by  structural  anisotropy,  which  only 
affected  the  constants. 

These  conclusions  are  In  general  agreement  with  those  reported 
by  Kopperman  et  al  (1982).  For  the  purposes  of  comparison,  results 
obtained  In  this  series  of  tests  and  those  reported  by  Kopperman  et  al 
(1982)  are  summarized  In  Table  4.2  and  shown  together  In  Fig.  4.3.  The 
values  of  slopes  obtained  In  the  present  series  of  tests  are  less  than 
those  reported  by  Kopperman  et  al  (1982).  However,  the  trend  Is  the  same 
In  both  cases  as  shown  In  Fig.  4.3,  namely  velocities  In  the  horizontal 
directions  are  higher  than  velocities  In  the  vertical  direction.  Illus¬ 
trating  the  effect  of  structural  anisotropy.  In  the  present  series  of 
tests,  velocities  In  the  horizontal  directions  (NS  and  EW)  are  almost  the 


TABLE  4.2 


Comparison  of  Constants  and  Slopes  for  Equation  4.2* 
with  Those  Reported  by  Kopperman  et  al  (1982) 


DIRECTION 

PRESENT  STUDY 

KOPPERMAN  et  al 

PRESENT  STUDY 

NS 

308 

250 

0.22 

EW 

280 

210 

0.23 

TB 

371 

271 

0.17 

.23 

.24 

.20 


E 


Confining  Pressure,  oq  (psi) 

Fig.  4.3  -  Comparison  of  P-wave  Velocities  with  Results 
from  Kopperman  et  al  (1982) 


same.  This  clearly  shows  the  Improvements  of  homogeneity  In  the  lateral 
directions  obtained  by  using  the  new  ralner  to  construct  the  sample  as 
discussed  In  Section  2.3. 


STRESS-STRAIN  BEHAVIOR  AND  CROSS-ANISOTROPIC  MODEL 


5.1  INTRODUCTION 

A  study  of  an  appropriate  constitutive  model  for  the  sand 
specimen  was  performed  In  light  of  the  stress-strain  data  provided  by 
the  stress  cells  and  strain  sensors  which  were  Installed  In  the  sand 
specimen  during  construction.  Since  the  sand  specimen  Is  Inherently 
anisotropic  (structural  anisotropy),  models  other  than  that  for  a 
homogeneous,  Isotropic  and  linear  elastic  material,  which  requires  two 
elastic  constants  (e.g.  Young's  modulus,  E,  and  Poisson's  ratio,  v), 
were  examined.  Based  on  the  results  of  P-wave  velocities  under  Isotrop¬ 
ic  confinement  (Section  4.3),  the  sand  specimen  can  be  idealized  as  a 
cross-anisotropic  material  (also  known  as  transversely  or  hexagonal ly 
anisotropic).  A  cross-anisotropic  material  Is  one  with  a  plane  of  sym¬ 
metry  In  which  all  stress-strain  relations  are  Isotropic:  l.e.,  the 
elastic  constants  are  the  same  no  matter  how  the  axes  are  chosen  within 
the  given  plane.  However,  the  elastic  constants  for  stresses  and 
strains  outside  this  plane  are  different.  Due  to  the  manner  In  which 
the  sand  was  rained  Into  the  trlaxlal  device  (Section  3.3),  the  hori¬ 
zontal  plane  turned  out  to  be  the  plane  of  symmetry.  This  Is  also  true 
for  many  naturally  deposited  soils  and  sedimentary  rocks. 


The  performances  of  stress  cells  and  strain  sensors  are 
examined  herein  In  relation  to  the  cross-anisotropic  model.  The 
stress-strain  data  and  wave  velocities  are  then  used  to  evaluate  the 
constants  for  the  cross-anisotropic  model. 


5.2  NOTATION  FOR  STRESSES,  STRAINS  AND  WAVE  VELOCITIES 


Conventional  notation  used  In  mechanics  Is  adopted  here; 
that  Is,  tension  is  considered  positive.  This  convention  Is  used 
because  compressive  strains  obtained  from  readings  of  the  strain  sen¬ 
sors  embedded  In  the  sand  specimen  gave  negative  values.  A  schematic 
drawing  of  stresses  acting  In  the  positive  sense  Is  shown  In  Fig.  5.1. 
It  should  be  noted  that  only  stresses  acting  on  the  exposed  faces  of 
the  element  are  shown  while  equal  complementary  stresses  acting  in 
opposite  directions  on  the  hidden  faces  are  omitted  for  simplicity.  The 
notation  used  for  the  stress  Is  for  the  first  subscript  to  denote  the 
plane  and  the  second  subscript  to  denote  the  line  of  action.  Therefore, 
normal  stresses  have  same  letters  for  the  two  subscripts  while  shear 
stresses  have  different  letters. 


Tensor  notation  Is  also  used.  Expressions  for  strains  are 
listed  In  the  following  equations  for  easy  reference: 


v.v.v.* •v-vv‘..'£. 

••  -■  .v.v  iv" 


exx  s  1,^3u/3x  +  3u/3x)  s  3u/3x  (5.1) 

eyy  =  *^3v^3y  +  3v/3y)  s  3v^3y  (5.2) 

E 22  =  *s(aw/3Z  +  3w/3z)  -  3W/3Z  (5.3) 

exy  =  eyx  =  ^(3v/3x  +  au/3y)  (5.4) 

eyz  =  e2y  =  >,(3W/3y  +  3V/3Z)  (5.5) 

ezx  =  exz  =  1*^3U/3Z  +  aw/3x)  (5.6) 


The  letters  u,  v  and  w  represent  displacements  In  the  x-,  y-  and 

z-dlrectlons,  respectively.  It  should  be  noted  that  shear  strains 

represented  by  exy»eyzancl  ezx  (Eqs.  5.4,  5.5  and  5.6)  are  exactly 

one-half  the  values  of  the  shear  strains  represented  by  y  ,  y  and 

xy  yz. 

yzx»  respectively. 

For  most  numerical  purposes,  stresses  and  strains  discussed 
above  are  written  as  vectors: 


to,T  ■  [oxx  °yy  °zz  °xy  V  “zx3 

(5.7) 

^  "  ^xx  eyy  Gzz  exy  eyz  Gzx^ 

(5.8) 

The  stress-strain  (constitutive)  laws  relating  the  two  are  usually 
expressed  In  matrix  form: 


where  [C]  and  [D]  are  6x6  matrices. 

Velocities  of  three  compression  waves  and  six  shear  waves 
propagating  and  polarized  along  principal  stress  directions  were  meas¬ 
ured  In  each  test.  These  nine  wave  velocities  are  denoted  by  the  sym¬ 
bols  shown  In  Fig.  5.2.  The  notation  of  velocity  used  here  Is  for  the 
first  subscript  to  denote  the  direction  of  wave  propagation,  while  the 
second  subscript  denotes  the  direction  of  polarization  (particle 
motion).  Since  the  direction  of  propagation  for  P-waves  Is  the  same  as 
the  direction  of  particle  motion,  P-wave  velocities  have  same  letters 
for  the  subscripts.  On  the  other  hand,  the  direction  of  propatatlon  for 
S-waves  Is  different  from  the  direction  of  particle  motion;  therefore, 
S-wave  velocities  have  different  letters  for  the  first  and  second  sub¬ 
scripts.  It  should  be  noted  that  the  xy-plane  Is  the  horizontal  plane 
which  Is  the  plane  of  symmetry  In  the  cross-anisotropic  model.  The  x-, 
y-  and  z-  axes  In  Figs.  5.1  and  5.2  are  the  principal  axes  In  the 
north-south  (NS),  east-west  (EW)  and  vertical  (TB)  directions,  respec¬ 
tively. 

The  above  notation  for  stresses,  strains  and  wave  velocities 
are  used  In  all  subsequent  analyses. 


5.3  PERFORMANCE  OF  STRESS  CELLS 


Three  total  stress  cells  (Section  2.5.5)  were  used  to  moni¬ 
tor  the  stresses  Induced  within  the  sand  mass  by  the  applied  confining 
pressures.  To  Investigate  any  effect  of  side  shear  along  the  walls  of 
the  trlaxlal  device  and  hence  to  check  the  uniformity  of  stresses  with¬ 
in  the  sand  mass,  two  of  the  three  stress  cells  (SS-1  and  SS-2  In  Fig. 
3.1)  were  located  such  that  they  both  measured  stresses  In  the 
x-dlrectlon.  The  third  stress  cell  (SS-3  In  Fig.  3.1)  was  used  to  meas¬ 
ure  stress  in  the  z-dl recti  on. 


Since  confining  pressures  applied  by  the  air-water  pressur¬ 
izing  system  (Figs.  1.6  and  1.7)  remained  stable  within  acceptable  lim¬ 
its  (±0.2  psl  (±1.4  kPa)),  stresses  Indicated  by  the  stress  cells 
before  the  top  membrane  leaked  were  plotted  against  the  applied  confin¬ 
ing  stresses  as  shown  in  Figs.  5.3  to  5.5.  Raw  readings  from  the  stress 
cells  together  with  calibration  curves  shown  in  Figs.  2.15  to  2.17  were 
used  to  obtain  stresses  from  the  stress  cells.  As  shown  In  Figs.  5.3  to 
5.5,  stresses  Indicated  by  the  stress  cells  were  widely  scattered  with 
all  but  four  readings  less  than  the  applied  pressure.  The  mean  values 
and  standard  deviations  were  calculated  for  stresses  obtained  under 
the  same  confining  stress.  Linear  regression  was  performed  on  the  mean 


the  Trlaxlal  Device 


Applied  Pressure  (pal) 

Cowpariaon  of  Principal  Straaa  Meaaured  by  Streaa  Call  SS-2  with  Applied  Streas  In 
the  Trlaxlal  Device 
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values  and  the  resulting  best-fit  straight  lines  are  shown  In  Figs.  5.3 
and  5.5. 


As  can  be  seen  In  the  figures,  readings  from  stress  cell  SS-2 
are  less  scattered  than  readings  from  the  other  two  stress  cells.  This 
Is  reasonable  since  cell  SS-2  Is  closer  to  the  pressurizing  membrane 
than  the  other  two  stress  cells  and  hence  less  Influenced  by  any  side 
shear  along  the  walls  of  the  trlaxlal  device.  However,  stresses 
obtained  from  cell  SS-1  are  larger  than  stresses  obtained  from  stress 
cell  SS-2.  This  Is  not  expected  since  side  shear  would  reduce  the 
transmitted  normal  stresses,  and  cell  SS-1  Is  located  farther  away  from 
the  pressurizing  membrane  than  cell  SS-2. 


Since  stresses  obtained  from  the  stress  cells  were  so  scat¬ 
tered,  It  was  concluded  that  further  effort  In  analyzing  the  stress 
data  was  unworthy.  Stress  readings  from  the  pressure  gauge  In  the 
pressurizing  system  were,  hence,  taken  to  be  representative  of  stress¬ 
es  In  the  sand  specimen  and  were  used  in  all  subsequent  analyses 
Involving  stresses. 


5.4  PERFORMANCE  OF  STRAIN  SENSORS 


Eight  pairs  of  strain  sensors  (Section  2.5.6)  were  Installed 
In  the  sand  specimen  In  an  attempt  to  measure  strains  Induced  by  the 
applied  confining  pressures.  The  locations  of  the  strain  sensors  are 
shown  in  Figs.  2.11  and  3.1.  Readings  of  null  amplitude  from  the  con¬ 
trol  unit  (Section  2.5.6)  for  each  pair  of  strain  sensors  were  recorded 
in  each  test.  These  amplitudes  together  with  the  calibration  curves 
shown  In  Fig.  2.18  and  Appendix  A  were  used  to  estimate  strains.  Howev¬ 
er,  strains  measured  by  strain  sensors  SN-1,  SN-2  and  SN-4  (Fig.  3.1) 
were  not  used  in  the  analyses  for  stress-strain  behavior  since  they  did 
not  measure  strains  in  the  principal  directions. 

It  was  soon  found  in  the  isotropic  series  of  tests  that  the 
strain  sensors  were  not  sensitive  enough  to  measure  strains  accurately 
since:  1.  null  amplitudes  changed  very  slightly  or  even  remained 
unchanged  under  increasing  or  decreasing  stress,  and  2.  for  some  cases, 
amplitude  readings  indicated  tensile  strains  under  increasing  isotrop¬ 
ic  confining  stress.  These  are  shown  in  Fig.  5.6.  Therefore,  amplitude 
readings  for  strains  were  not  measured  in  most  of  the  tests  under  biax¬ 
ial  confinement  before  the  top  membrane  leaked.  However,  amplitude 
readings  were  again  measured  after  a  new  top  membrane  was  Installed 
hoping  that  some  useful  strain  data  might  be  obtained.  Strain  readings 
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were  used  to  calculate  Poisson's  ratios and  v,.u  where  the  sub- 

HV  VH 


scripts  H  and  V  denote  the  horizontal  and  vertical  directions,  respec¬ 
tively.  The  relation  between  horizontal  strains  and  vertical  loadings 
or  between  vertical  strains  and  horizontal  loadings  is  controlled  by 


v  or  v.#u  •  However,  It  should  be  noted  that  strain  data  were  ana- 
H  V  VH 


lyzed  with  biased  judgement  since  less  than  a  quarter  of  the  data  was 
used  In  analyses.  Therefore,  calculations  using  strain  data  only  gave 
rough  approximations  of  the  constants. 


5.5  CROSS-ANISOTROPY 


The  constitutive  model  for  the  sand  specimen  can  be 
expressed  by  Eq.  5.9  or  5.10.  Due  to  cross-anisotropy,  there  are  only 
five  Independent  constants  to  be  evaluated  In  the  stress-strain  matrix 
[C]  of  Eq.  5.9  or  strain-stress  matrix  [D]  of  Eq.  5.10  : 
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Shear  wave  velocities,  In  addition  to  normal  stress-strain 
data,  had  to  be  used  to  evaluate  the  five  Independent  constants  In  the 
matrix  [D]  since  shear  stresses  and  strains  were  not  measured.  Once 
the  terms  In  matrix  [D]  were  evaluated,  matrix  [C]  could  be  obtained  by 
Inverting  matrix  [D].  On  the  other  hand,  the  nine  wave  velocities 
(Section  5.2)  measured  along  principal  directions  could  be  used  to 
evaluate  four  of  the  five  constants  in  matrix  [C]  and  to  estimate  the 
fifth  constant.  These  two  approaches  are  discussed  In  the  following 


sections. 


5.5.1  Evaluation  of  Constants  by  Stress-Strain  Data  and  Wave 
Velocities 

For  the  cross-anisotropic  model  with  the  xy-plane  (horizon¬ 
tal  plane)  as  the  plane  of  symmetry,  the  terms  In  the  strain-stress 
matrix  [Dl  In  Eq.  5.12  as  suggested  by  Desal  and  Christian  (1977)  are: 


The  symbols  E^  and  Ey  represent  Young's  moduli  In  the  hori¬ 
zontal  plane  and  vertical  direction,  respectively.  Poisson's  ratio 
relating  the  strains  In  one  horizontal  axis  to  strains  In  the  other 
horizontal  axis  Isv^  while  VyH  and  v^are  defined  previously  In  Sec¬ 
tion  5.4  and  are  not  Independent  of  each  other  because  of  symmetry  In 
the  [0]  matrix.  The  Gy^term  relates  shear  stresses  out  of  the  horizon¬ 
tal  plane  to  shear  strains  out  of  the  horizontal  plane. 

Using  the  terms  In  Eq.  5.13  for  the  [0]  matrix  and  expanding 
the  matrices  In  Eq.  5.10  for  normal  strains  gives: 


It  should  be  noted  that  for  biaxial  confinements  with  same 


horizontal  stresses,  Eq.  5.14  Is  the  same  as  Eq.  5.15  and  strains  In 
the  x-dlrectlon  should  be  equal  to  the  strains  In  the  y-dlrectlon. 
However,  strain  readings  from  the  strain  sensors  gave  values  that  were 
far  off  from  each  other,  up  to  an  order  of  magnitude  different  showing 
the  poor  quality  of  strain  data.  Some  usable  strain  data  under  appro** 
prlate  biaxial  conditions  gave  average  values  of  0.35  and  0.17  for  v^y 
and  VyH  respectively.  However,  It  should  be  emphasized  again  that 
these  were  only  rough  estimations  for  the  parameters  since  the  strain 
data  were  not  satisfactory  as  discussed  In  Section  5.4.  The  value  of 
Poisson's  ratio  v^was  calculated  using  the  following  equation: 
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With  data  of  S-wave  velocities  (Lee  et  al,  1984)  and  P-wave  velocities 
for  waves  propagating  along  horizontal  axes  under  Isotropic  confine¬ 
ment,  the  average  value  of  V^was  found  to  be  0.27.  With  known  values 
for  the  Poisson's  ratios,  Eqs.  5.14  to  5.16  were  solved  for  E^  and  Ey 
using  the  available  stress-strain  data.  Due  to  the  poor  quality  of  the 
data,  values  of  Eu  ranged  from  9,300  psl  (64  MPa)  to  70,400  psl  (485 
MPa)  with  an  average  of  21,100  psl  (145  MPa)  while  values  of  Ey  ranged 
from  7,700  psl  (53  MPa)  to  9,200  psl  (63  MPa)  with  an  average  of  8,300 
psl  (57  MPa).  Shear  wave  velocities  provided  from  Lee  et  al  (1984) 
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together  with  equation  suggested  by  Bachman  (1983)  were  used  to  evalu¬ 
ate  the  last  parameter  .  The  equation  used  was: 

GVH  =  [e(pVzx2)  +  E(pvzy2)  +  z(pVyz2)  +  E(pVxz2)]/n  (5.18) 

where  n  Is  the  number  of  measurements  and  pis  the  mass  density  of  the 
sample.  The  average  value  of  Gy^ calculated  from  shear  wave  velocities 
under  Isotropic  confinements  at  20  psl  (138  kPa)  was  17,700  psl  (122 
MPa). 


As  suggested  by  Desal  and  Christian  (1977),  the  terms  In  the 
cross-anisotropic  [C]  matrix  found  by  Inverting  the  [D]  matrix  are: 

Cn  -  C22  =  A(|*-  vJH)  =  ^ 

C33  =  =  Mv 

c“4  ■  *  2shh 

C55  =  C66  *  2GyH  (5.19) 

eh 

C12  =  C21  =  A{|^vhh  +  vJH)  =  cn  -  C44 
cl3  *  C31  =  C23  =  C32  =  AvyH(l+vHn) 
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Since  the  values  of  the  parameters  In  the  [D]  matrix  were 
already  calculated,  the  terms  In  the  [C]  matrix  had  the  following  val¬ 
ues: 
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It  should  be  noted  that  the  ratio  of  to  My  (2.70)  Is  too 
large  when  compared  to  the  square  of  the  ratio  of  P-wave  velocity  In 
the  horizontal  direction  to  P-wave  velocity  In  the  vertical  direction, 
the  value  of  which  as  approximated  by  Eq.  4.1  Is  1.37. 


5.5.2  Estimation  of  Constants  by  Wave  Velocities 

I 

With  measurements  of  the  nine  wave  velocities  under  isotrop¬ 


ic  confinement,  four  of  the  five  Independent  parameters  In  the  [C] 
matrix  can  be  evaluated.  The  equations  suggested  by  Bachman  (1983)  to 
evaluate  these  parameters  are: 


Cll  =  C22  a  ^  s  [E(pVxx2)  +  S(pVyy2)]/ n 


C33  =  My  =  z(pVzz2)/n 
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The  expressions  for  the  remaining  terms  In  the  [C]  matrix  suggested  by 
Drnevlch  (1975)  are: 

C12  s  C21  *  %  -  2(5^  (5.21) 

C13  =  C23  =  C32  =  C31  -  «s(^  +  My)  -  2Gvh  (5.22) 

It  should  be  noted  that  Eq.  5.22  Is  only  an  approximation  since  meas¬ 
urement  of  P-wave  velocities  propagating  at  an  angle  to  the  principal 
axis  Is  needed  to  evaluate  the  fifth  constant.  Equation  for  calculat¬ 
ing  the  fifth  constant  from  oblique  P-waves  as  presented  by  Podlo 
(1968)  Is: 

c13  •  { CZp V^5 „  -  MC,,  +  C33  ♦  2C44)]2  -  [Cu  ~  ‘S2}'1  -  C44  (5.23) 

where  V  o  is  the  velocity  of  P-wave  propagating  at  45  degrees  to  the 
principal  directions. 
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With  S-wave  velocities  (Lee  et  al  ,  1984)  and  P-wave  veloci¬ 
ties  (Table  4.1)  under  Isotropic  confinement  of  20  psl  (138  kPa),  the 
expressions  In  Eq.  5.20  are  found  to  have  the  following  values: 

Cn  =  64100  psi  (442  MPa) 

Cs3  =  48000  psi  (331  MPa) 

°44  =  46500  psi  (320  MPa) 

C55  *  35400  psi  (244  MPa) 

From  the  above  results,  the  remaining  terms  In  the  [C]  matrix  calcu¬ 
lated  by  Eqs.  5.21  and  5.22  are: 

C12  =  17700  psi  (122  MPa) 

C13  »  20700  psi  (143  MPa) 

When  compared  to  values  of  the  same  terms  determined  In  Sec¬ 
tion  5.5.1,  the  above  terms  are  larger.  The  largest  difference 
occurred  In  the  term  C13  with  a  value  about  10  times  larger. 

5.6  SUMMARY 

The  stress-strain  behavior  of  the  sand  specimen  Is  studied. 

To  help  In  this  study,  notation  for  stresses,  strains  and  wave  veloci¬ 
ties  are  first  defined  In  Section  5.2.  The  performance  of  stress  cells 
and  strain  sensors  are  then  evaluated.  Unfortunately,  both  stress 


cells  and  strain  sensors  are  found  to  have  provided  data  which  were 
highly  scattered  and  generally  unsatisfactory. 

The  cross-anisotropic  model  for  the  sand  specimen  is  studied 
analytically  In  Section  5.5.  With  available  readings  of  stresses, 
strains  and  wave  velocities,  complete  evaluation  of  all  constants  In 
the  cross-anisotropic  matrix  [C]  relating  stress  and  strain  can  be 
accomplished.  However,  It  should  be  noted  that  data  of  stresses, 
strains  and  wave  velocities  were  obtained  from  different  stress 
states.  It  was  also  found  that  a  consistent  set  of  wave  velocity  (both 
P-wave  and  S-wave)  data  could  be  used  to  evaluate  four  of  the  five 
Independent  constants  In  the  [C]  matrix  and  to  approximate  the  fifth. 
If  oblique  P-waves  were  monitored  and  velocities  were  measured  In  the 
present  series  of  tests,  the  fifth  constant  could  have  been  evaluated 


by  these  velocities. 


CHAPTER  SIX 


BIAXIAL  CONFINEMENT 


6.1  INTRODUCTION 

After  studying  the  effect  of  Isotropic  confinement  on  P-wave 
velocity,  an  exhaustive  set  of  tests  with  biaxial  confinement  was  per¬ 
formed.  Conditions  of  biaxial  confinement  were  obtained  by  keeping  the 
intermediate  principal  effective  stress  equal  to  the  minor  principal 
effective  stress  (01>52=03)  or  ^eeP^n9  the  Intermediate  principal 
effective  stress  equal  to  the  major  principal  effective  stress 
(ct1=ct2>ct3). 


The  complete  set  of  tests  was  composed  of  two  series  of  biax¬ 
ial  confinement  tests:  the  first  series  consisted  of  tests  with  confin¬ 
ing  stress  varying  In  only  one  principal  direction  while  the  second 
series  consisted  of  tests  with  confining  stress  varying  in  two  princi¬ 
pal  directions.  Additionally,  the  first  series  of  tests  contained  two 
subsets:  the  first  subset  was  composed  of  tests  with  the  Intermediate 
principal  stress  always  equal  to  the  minor  principal  stress  while  the 
second  subset  was  composed  of  tests  where  the  Intermediate  principal 
stress  ranged  In  value  from  the  minor  to  the  major  principal  stress. 
The  loading  conditions  for  this  first  series  of  tests  are  shown  In 


Figs.  6.1  and  6.2  for  subsets  one  and  two,  respectively.  The  loading 
conditions  for  the  second  test  series  are  shown  In  Fig.  6.3.  Tests  of 
unloading  sequences  for  the  conditions  shown  in  Figs.  6.1  to  6.3  were 
also  performed  so  that  complete  data  were  available  for  the  study  of 
stress  history  and  possible  structural  anisotropy  effects. 

Tests  with  biaxial  confinement  started  with  an  Isotropic 
state  of  stress  of  15  psl  (103  kPa).  The  stress  in  the  vertical  direc¬ 
tion  (z-directlon)  was  then  Increased  from  15  psl  (103  kPa)  to  20,  30 
and  40  psi  (138,  207  and  276  kPa)  as  shown  In  Fig.  6.1(a).  With  the 
stress  of  40  psl  (276  kPa)  being  held  constant  In  the  z-direction, 
stresses  In  the  x-  and  y-  directions  were  then  increased  from  15  to  40 
psi  (103  to  276  kPa)  in  the  same  Increments  as  before.  This  loading 
condition  Is  shown  In  Fig.  6.3(a).  Then  the  horizontal  stresses  were 
unloaded  from  40  psi  (276  kPa)  to  15  psi  (103  kPa)  in  the  reverse 
sequence  with  the  vertical  stress  staying  constant  at  40  psl  (276  kPa). 
With  the  horizontal  stresses  being  held  constant  at  15  psi  (103  kPa), 
the  vertical  stress  was  unloaded  from  40  psi  (276  kPa)  to  15  psl  (103 
kPa)  In  the  same  reverse  sequence.  This  represented  major  principal 
stress  In  the  z-direction.  Then,  complete  series  of  tests  was  redone 
with  5i  In  the  x-dlrection.  The  loading  conditions  are  shown  in  Figs. 
6.1(b)  and  6.3(b).  Then,  complete  series  of  tests  was  redone  with  5i  in 


(a)  First  Orientation  of  Principal  Stresses 


a  -  15  psi 


(b)  Second  Orientation  of  Principal  Stresses 


a 


a  =  15  psi 


15  -  40  psi 


(c)  Third  Orientation  of  Principal  Stresses 


Fig.  6.1  -  Initial  Subset  of  First  Series  of 
Biaxial  Confinement  Tests  with 
Variation  of  Stress  in  One  Principal 
Direction 


(a)  First  Orientation  of  Principal  Stresses 


o  =  20  psi 


(b)  Second  Orientation  of  Principal  Stresses 

Fig.  6.2  -  Second  Subset  of  First  Series  of 
Biaxial  Confinement  Tests  with 
Variation  of  Stress  in  One  Principal 
Direction 


the  y-dlrectlon  with  loading  conditions  shown  In  Figs.  6.1(c)  and 
6.3(c). 

In  addition  to  the  tests  described  above,  tests  with  loading 
conditions  shown  In  Fig.  6.2  were  performed  to  study  the  effect  of 
Intermediate  principal  effective  stress  on  P-wave  velocity.  As  shown 
In  Fig.  6.2(a),  stress  In  the  z-dlrectlon  ranged  In  value  from  the 
minor  to  the  major  principal  stress  while  stresses  In  the  x-  and  y- 
dlrectlons  remained  constant  at  20  psl  (138  kPa)  which  represented  the 
Intermediate  principal  stress.  It  should  be  noted  that  the  Intermedi¬ 
ate  principal  stress  was  also  the  major  principal  stress  when  the 
stress  In  the  z-dlrectlon  was  the  minor  principal  stress,  and  vice  ver¬ 
sa.  Then,  complete  series  of  tests  was  redone  with  stress  variation  in 
the  x-dlrectlon  as  shown  In  Fig.  6.2(b).  Tests  with  stress  variation 
in  the  y-dlrectlon  are  not  performed  since  the  sample  Is  believed  to  be 
cross-anisotropic  (Section  5.1)  and  enough  data  were  believed  to  have 
been  obtained  for  studying  the  effect  of  Intermediate  principal  stress 
on  P-wave  velocities. 


The  exact  sequences  of  the  complete  series  of  tests  with 
biaxial  confinement  are  summarized  In  Table  3.2  and  shown  schemat¬ 


ically  In  Figs.  3.6  and  3.7. 


6.2  VARIATION  OF  STRESS  IN  ONE  PRINCIPAL  DIRECTION 


Compression  wave  velocities  were  measured  along  each  princi¬ 
pal  stress  axis  under  each  biaxial  state.  The  same  procedure  discussed 
In  Section  3.5  was  used  to  analyze  all  data.  The  biaxial  conditions 
shown  In  Figs.  6.1  and  6.2  are  referred  to  as  the  first  series  of  tests 
with  variation  of  stress  In  one  principal  direction.  Results  from  the 
first  subset  of  these  tests  (Fig.  6.1)  were  used  to  study  the  effects 
of  stress  history,  structural  anisotropy  and  biaxial  confinement  on 
P-wave  velocities.  Results  from  the  second  subset  of  the  first  test 
series  (Fig.  6.2)  were  used  to  study  the  effect  of  Intermediate  princi¬ 
pal  effective  stress  on  P-wave  velocities. 

6.2.1  Effect  of  Stress  History 

Compression  wave  velocities  obtained  under  the  same  loading 
and  unloading  conditions  for  each  biaxial  state  of  stress  were  used  to 
study  the  effect  of  stress  history  on  P-wave  velocities.  One  compli¬ 
cation  was  that  the  unloading  conditions  were  performed  sometime  after 
each  of  the  corresponding  loading  conditions  shown  In  Fig.  6.1.  This 
time  delay  was,  however,  assumed  Inconsequential. 


Compression  wave  velocities  propagating  along  the  axis  of 
Increasing  stress  are  plotted  against  the  corresponding  stress  In  Fig. 
6.4  for  both  loading  and  unloading  conditions.  The  notation  for  veloc¬ 
ities  discussed  in  Section  5.2  Is  used.  The  notation  for  stresses  used 
by  Kopperman  et  al  (1982)  Is  also  used  In  which: 

a  =  effective  stress  In  the  direction  of  wave  propagation, 

d 

0(3  =  effective  stress  In  the  direction  of  particle  motion,  and 
5C  =  effective  stress  perpendicular  to  agand  5^ 

For  P-waves,  afl  and  5^  are  In  the  same  direction.  (For  shear  waves, 
5a  and  5^  are  perpendicular.) 

As  shown  In  Fig.  6.4,  P-wave  velocities  in  the  unloading  con¬ 
dition  are  slightly  higher  than  those  In  the  loading  condition. The  dif¬ 
ferences  are  small  with  the  largest  discrepancy  of  about  8  percent 
occurring  at  20  psl  (138  kPa)  In  the  y-di recti  on.  The  mean  value  of  the 
percentage  differences  for  the  rest  of  the  data  In  all  three  directions 
Is  2  percent.  Hence,  the  effect  of  stress  history  Is  negligible.  As  a 
result,  average  values  of  the  data  from  loading  and  unloading  are  used 
for  all  subsequent  analyses. 


O  loading 
A  unloading 
m  =  slope 
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Fig.  6.4  -  Effect  of  Stress  History  on  P-wave  Velocity 
Under  Biaxial  Confinement  with  Variation  of 
Stress  in  One  Direction 


6.2.2  Effect  of  Structural  Anisotropy 


If  the  sand  specimen  was  perfectly  uniform  and  Isotropic, 
P-wave  velocities  shown  In  Fig.  6.4  along  the  axis  of  Increasing  stress 
would  have  been  the  same  for  each  of  the  three  different  orientations 
of  the  major  principal  stress.  Best-fit  straight  lines  for  the  average 
values  of  the  loading  and  unloading  data  for  each  of  these  three  biaxi¬ 
al  conditions  are  shown  together  In  Fig.  6.5.  The  three  straight  lines 
show  exactly  the  same  trend  as  those  for  Isotropic  conditions  (Fig. 
4.2).  Therefore,  the  effect  of  structural  anisotropy  Is  Independent  of 
these  stress  states  and  has  the  same  effect  on  P-wave  velocities  for 
both  Isotropic  and  biaxial  conditions. 

As  shown  In  Fig.  6.5,  the  cross-anisotropic  model  for  the 
sand  specimen  Is  again  confirmed  since  P-wave  velocities  In  the  hori¬ 
zontal  directions  (NS  and  EW)  are  almost  the  same,  with  a  maximum  per¬ 
centage  difference  of  less  than  4  percent.  Due  to  cross  anisotropy, 
velocities  In  the  horizontal  directions  are  about  17  percent  higher 
than  the  velocity  In  the  vertical  direction  for  the  same  state  of 
stress. 


Confining  Stress  in  the  Direction  of  Wave  Propagation 
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Fig.  6.5  -  Effect  of  Structural  Anisotropy  on  P-wave 
Velocities  Under  Biaxial  Confinement  with 
Variation  of  Stress  in  One  Direction 
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6.2.3  Effect  of  Biaxial  Confinement 

To  study  the  effect  of  biaxial  confinement  onVp,  velocities 
of  waves  propagating  along  each  principal  axis  for  conditions  shown  In 
Fig.  6.1  are  plotted  against  the  Increasing  confining  stress  In  Figs. 
6.6,  6.7  and  6.8  respectively.  Compression  wave  velocities  under  Iso¬ 
tropic  conditions  are  also  shown  In  the  figures  for  comparison  pur¬ 
poses. 

As  shown  In  Figs.  6.6(a),  6.7(a)  and  6.8(a),  compression 

wave  velocity  In  the  direction  of  Increasing  confining  stress  under 

biaxial  confinement  varies  with  stress  In  essentially  the  same  manner 

as  under  isotropic  confinement.  These  results  suggest  that  the 

relationship  between  wave  velocity  and  confining  stress  for  biaxial 

confinement  can  be  patterned  after  Eq.  4.2  for  Isotropic  confinement  If 

5  (stress  In  the  direction  of  wave  propagation  in  psf)  Is  substituted 
d 

for  dQ  (Isotropic  confining  stress  In  psf).  The  resulting  equation  Is: 

V„  -  (6.1) 

Values  of  constants  In  Eq.  6.1  were  obtained  by  fitting  a 
power  curve  to  each  set  of  data  using  least-squares  method  (Figs.  6.6a, 
6.7a  and  6.8a).  These  values  along  with  those  determined  for  isotropic 
confinement  are  summarized  In  Table  6.1.  Host  importantly,  since  the 


Biaxial,  m  =  0.20 


Isotropic,  m 


Biaxial,  m 


Confining  Stress  in  the  y-direction,  psi 

6.8  -  Variation  in  P-wave  Velocities  Under  Biaxial 
Conditions  for  Stress  Increasing  Along  y- 
direction 


equation  for  biaxial  confinement  only  relates  P-wave  velocity  to  the 
stress  in  the  direction  of  wave  propagation,  stresses  In  the  other  two 
principal  directions  perpendicular  to  the  direction  of  propagation 
have  very  little  effect  on  P-wave  velocity. 

As  shown  in  the  middle  and  bottom  curves  of  Figs.  6.6  to  6.8, 
P-wave  velocities  along  directions  of  constant  stress  are  almost  con¬ 
stant  with  the  largest  percentage  difference  less  than  6  percent.  This 
is  just  experimental  scatter  since  some  velocities  increase  and  some 
decrease.  The  results  are  reasonable  (for  velocities  to  remain  con¬ 
stant  in  the  direction  of  constant  stress)  If  P-wave  velocity  can  be 
expressed  by  Eq.  6.1. 

It  should  be  noted  that  at  15  psi  (103  kPa),  the  loading  con¬ 
dition  is  actually  an  Isotropic  confinement.  Therefore,  the  P-wave 
velocities  measured  under  this  loading  condition  should  be  the  same  as 
those  under  isotropic  confinement  and  the  two  straight  lines  should 
Intersect  each  other  at  15  psi  (103  kPa).  This  is  shown  In  Figs.  6.6  to 
6.8  with  small  experimental  scatter  and  once  again  the  effect  of  stress 
history  Is  confirmed  to  be  negligible. 

It  Is  also  Interesting  to  observe  that  complete  anisotropy 
of  the  system  resulted  from  the  coupling  of  structural  anisotropy  and 
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stress  anisotropy.  The  sand  specimen  Is  known  to  be  cross-anisotropic 
due  to  structural  anisotropy  (Section  4.3).  Under  biaxial  loading  with 
different  horizontal  stresses,  soil  properties  in  the  two  horizontal 
directions  will  also  be  different  due  to  stress-induced  anisotropy 
and,  therefore,  complete  anisotropy  resulted.  Complete  anisotropy 
under  biaxial  loading  is  shown  by  condensing  the  solid  curves  in  Figs. 
6.7(a)  to  6.7(c)  into  a  single  figure  which  is  shown  in  Fig.  6.9.  Com¬ 
pression  wave  velocities  in  each  direction  are  different  from  one 
another  (complete  anisotropy)  except  at  15  psi  (103  kPa)  where  the 
velocities  in  the  horizontal  directions  are  equal  for  the  isotropic 
confinement  (no  stress  anisotropy).  Therefore,  complete  anisotropy 
resulted  in  biaxial  conditions  shown  in  Figs.  6.1(b),  6.1(c),  6.2(b), 
6.3(b)  and  5.3(c). 

6.2.4  Effect  of  Intermediate  Principal  Effective  Stress 

The  results  presented  so  far  are  for  biaxial  conditions  with 
the  intermediate  principal  effective  stress  always  equal  to  the  minor 
principal  effective  stress.  For  the  biaxial  conditions  shown  in  Fig. 
6.2,  the  Intermediate  principal  stress  is  equal  to  the  major  principal 
stress  at  the  start  of  testing  and  then  changes  to  be  equal  to  the 
minor  principal  stress  in  the  later  stages.  Tests  with  these  biaxial 
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Confining  Stress  in  the  x-direction,  psi 


Fig.  6.9  -  Effect  of  Complete  Anisotropy  on  P-wave 
Velocities  Under  Biaxial  Condition  with 
Different  Horizontal  Confining  Stresses 
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conditions  were  performed  to  Investigate  If  there  Is  any  effect  of 
Intermediate  principal  stress  on  P-wave  velocity. 


Results  of  the  biaxial  conditions  shown  In  Figs.  6.2(a)  and 
6.2(b)  are  shown  In  Figs.  6.10  and  6.11,  respectively.  Comparable 
biaxial  conditions  with  increasing  confining  stress  In  the  same 
directions  are  also  shown  In  Figs.  6.10  and  6.11.  Since  the  P-wave 
velocities  are  very  close  to  each  other  for  both  stress  conditions, 
there  Is  no  discernible  effect  of  Intermediate  principal  stress  on  the 
P-wave  velocities  and  Eq.  6.1  Is  applicable  for  both  cases.  Com¬ 
pression  waves  propagating  along  directions  of  constant  stress  have 
almost  constant  values  for  P-wave  velocities  and  the  results  appeared 
as  horizontal  straight  lines  In  Figs.  6.10  and  6.11. 

6.3  VARIATION  OF  STRESS  IN  TWO  PRINCIPAL  DIRECTIONS 

For  the  conditions  shown  In  Fig.  6.3,  two  of  the  three  prin¬ 
cipal  stresses  were  changed  during  each  test.  Velocities  of  com¬ 
pression  waves  propagating  along  axes  with  increasing  confining  stress 
are  plotted  against  the  confining  stress  In  Fig.  6.12.  Compression 
wave  velocities  obtained  under  Isotropic  confinement  are  also  shown 
for  comparison  purposes.  All  of  the  resulting  straight  lines  are  In 
close  agreement  to  each  other  showing  same  trend  In  the  variation  of  Vp 
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Confining  Stress  in  the  x-direction,  psi 

Fig.  6.11  -  Effect  of  Intermediate  Principal  Stress  on 
V_  with  Stress  Increasing  Along  the  x- 
dlrection 
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with  confining  stress.  The  largest  difference  between  P-wave  veloci¬ 
ties  under  Isotropic  confinement  and  biaxial  confinement  Is  shown  in 
Fig.  6.12(c)  at  15  psl  (103  kPa)  at  which  point  the  difference  Is  12 
percent.  The  average  of  the  percentage  differences  In  P-wave  veloci¬ 
ties  shown  In  Fig.  6.12  Is  4  percent.  The  scatter  Is  small  and  once 
again  Eq.  6.1  Is  applicable  for  P-wave  velocities  propagating  along 
axes  of  Increasing  confining  stress.  Values  of  the  constants  in  Eq. 
6.1  derived  from  the  results  presented  In  Fig.  6.12  are  summarized  in 
Table  6.2. 


Velocities  of  P-waves  propagating  along  axes  of  constant 
confining  stress  are  shown  In  Fig.  6.13.  The  P-wave  velocities  in  each 
of  the  principal  axes  with  constant  confining  stress  remained  more  or 
less  constant.  The  resulting  straight  lines  are  almost  horizontal  with 
small  values  for  the  slopes.  Structural  anisotropy  Is  also  depicted  In 
Fig.  6.13  since  the  velocities  In  the  horizontal  directions  are  almost 
the  same  and  are  about  17  percent  higher  than  the  velocities  In  the 
vertical  direction. 

6.4  SUMMARY 

From  the  results  of  P-wave  velocities  measured  under 
repeated  biaxial  conditions,  It  was  shown  that  the  effect  of  stress 


TABLE  6.2 


Comparison  of  Constants  and  Slopes  for  Equations 
Relating  Vp  to  5a  for  Isotropic  Confinement* 
and  Vp  to  3a  for  Biaxial  Confinement** 
with  Variation  of  Stress  in  Two  Principal  Directions*** 


DIRECTION 

STRESS  STATE 

Cl 

m 

ISOTROPIC 

308 

0.22 

X 

BIAXIAL  2+ 

431 

0.18 

BIAXIAL  3 

342 

0.20 

ISOTROPIC 

280 

0.23 

Y 

BIAXIAL  1 

366 

0.21 

BIAXIAL  3 

256 

0.24 

ISOTROPIC 

371 

0.17 

Z 

BIAXIAL  1 

227 

0.22 

BIAXIAL  2 

443 

0.16 

*V  ■  C,om  (Isotropic  Confinement,  Eq.  4.2) 
p  i  o 

**V  =  C  am  (Biaxial  Confinement,  Eq.  6.1) 

p  JL  3. 

***Biaxial  conditions  are  shown  in  Fig.  6.3 


+  See  insert  in  Fig.  6.12 
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history  on  P-wave  velocities  of  this  sand  was  negligible.  Therefore, 
averages  of  P-wave  velocities  measured  under  loading  and  unloading 
conditions  were  used  for  analysis  purposes,  and  the  small  effect  of 
stress  history  was  neglected. 

For  biaxial  conditions  with  equal  horizontal  confining 
stresses,  the  sand  specimen  behaved  as  a  cross-anisotropic  material 
with  wave  propagation  velocities  in  the  horizontal  directions  equal 
and  different  from  Vp  in  the  vertical  direction.  However,  when  the 
two  horizontal  confining  stresses  are  different,  the  resulting  stress 
anisotropy  couples  with  structural  anisotropy  to  produce  a  complete 
anisotropic  material.  Compression  wave  velocities  In  a  complete  aniso¬ 
tropic  material  have  different  values  in  each  of  the  three  principal 
axes.  However,  the  same  equation  (Eq.  6.1)  is  applicable  to  velocities 
of  P-waves  propagating  along  directions  of  Increasing  stress  In  the 
complete  anisotropic  state. 

Curves  for  average  values  of  velocities  of  P-waves  propagat¬ 
ing  along  each  principal  axis  with  Increasing  confining  stress  from  all 
the  appropriate  biaxial  conditions  shown  In  Figs.  6.1  to  6.3  are  shown 
In  Fig.  6.14.  The  three  straight  lines  are  almost  Identical  to  those 
In  Fig.  4.2  for  the  Isotropic  condition  except  that  P-wave  velocities 
In  the  vertical  direction  are  3  to  4  percent  less  than  those  under  Iso- 


wave  Velocity 


tropic  conditions.  The  values  of  constants  In  Eq.  6.1  for  the  straight 
lines  shown  In  Fig.  6.14  are  summarized  In  Table  6.3  together  with 
those  for  biaxial  loading  conditions  reported  by  Kopperman  et  al 
(1982). 
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CHAPTER  SEVEN 


TRIAXIAL  CONFINEMENT 


7.1  INTRODUCTION 


The  last  step  in  the  study  of  the  effect  of  confining  pres¬ 
sure  on  compression  wave  velocity  was  to  perform  tests  with  tri axial 
confinement.  Conditions  of  trlaxlal  confinement  were  obtained  when  all 
stresses  In  principal  directions  had  different  values,  l.e.  when  the 
major  principal  stress  was  greater  than  the  Intermediate  principal 
stress  which  was  In  turn  greater  than  the  minor  principal  stress 
(51>52  >^3 ) * 


The  complete  set  of  trlaxlal  confinement  tests  was  composed 
of  three  series  of  tests:  1.  the  first  series  consisted  of  tests  in 
which  confining  stress  was  varied  In  only  one  principal  direction,  2. 
the  second  series  consisted  of  tests  In  which  confining  stress  was  var¬ 
ied  In  two  principal  directions,  and  3.  the  third  series  consisted  of 
tests  In  which  confining  stress  was  varied  In  all  three  principal 
directions.  The  loading  conditions  for  each  series  are  shown  In  Fig. 
7.1.  The  exact  sequence  of  these  loading  conditions  Is  shown  In  Fig. 
3.7.  Unloading  tests  were  not  performed  since  results  from  previous 
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tests  (Sections  4.2  and  6.2.1)  had  shown  that  the  effect  of  stress  his¬ 
tory  on  P-wave  velocities  was  negligible. 

Compression  wave  velocities  were  measured  along  each  princi¬ 
pal  stress  axis  under  each  tri axial  state.  The  same  procedure  dis¬ 
cussed  In  Section  3.5  was  used  to  analyze  all  data.  Results  and 
conclusions  for  the  first,  second  and  third  series  of  triaxlal  confine¬ 
ment  tests  are  presented  In  Sections  7.2,  7.3  and  7.4,  respectively. 

7.2  VARIATION  OF  STRESS  IN  ONE  PRINCIPAL  DIRECTION 

As  shown  in  Fig.  7.1(a),  the  first  series  of  triaxlal  con¬ 
finement  tests  started  with  stresses  In  the  x-,  y-  and  z-dl recti ons 
having  values  of  15  psl  (103  kPa),  20  psi  (138  kPa)  and  40  psl  (276 
kPa),  respectively.  The  Intermediate  principal  stress  was  then 
increased  from  20  psl  (138  kPa)  to  35  psl  (241  kPa)  in  Increments  of  5 
psl  (34  kPa)  while  the  major  and  minor  principal  stresses  In  the  z-  and 
x-  directions,  respectively,  remained  unchanged.  Velocities  obtained 
with  these  stress  conditions  could  be  used,  therefore,  to  study  the 
effect  of  Intermediate  prirclpal  stress  on  P-wave  velocities  under 
triaxlal  confinement. 
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7.2.1  Effect  of  Intermediate  Principal  Stress 

Compression  wave  velocities  obtained  from  the  stress  condi¬ 
tions  shown  in  Fig.  7.1(a)  are  plotted  against  the  stress  in  the 
y-direction  (the  intermediate  principal  stress  in  this  case)  in  Fig. 
7.2.  As  shown  in  the  figure,  compression  waves  propagating  along  axes 
of  constant  stress  (x-  and  z-  axes)  travelled  at  almost  constant  values 
of  velocities  and,  hence,  were  independent  of  the  increasing  Interme¬ 
diate  principal  stress.  The  resulting  curves  therefore  appear  as  near¬ 
ly  horizontal  lines  in  Fig.  7.2.  On  the  other  hand,  P-waves  propagating 
along  the  axis  of  the  intermediate  principal  stress  exhibited  veloci¬ 
ties  which  increased  with  Increasing  stress  with  a  linear 
relationship  as  shown  in  Fig.  7.2.  As  a  result,  Eq.  6.1  (which  was  pat¬ 
terned  after  Eq.  4.2  for  isotropic  confinement)  for  biaxial  confine¬ 
ment  could  be  used  to  express  the  relationship  between  P-wave  velocity 
and  confining  stress  in  the  direction  of  wave  propagation  for  triaxial 
confinement. 

The  best-fit  straight  line  obtained  from  measurements  under 
isotropic  conditions  of  P-waves  propagating  In  the  y-dlrectlon  is  also 
shown  In  Fig.  7.2.  Since  the  two  straight  lines  are  almost  on  top  of 
each  other,  the  relationship  between  P-wave  velocity  and  confining 


Confining  Stress  in  the  yy-direction,  psi 

7.2  -  Variation  of  Compression  Wave  Velocities 
Under  Triaxial  Confinement  When  Only  One 
Principal  Stress  Varied 


stress  In  the  direction  of  wave  propagation  is  the  same  for  both  iso¬ 
tropic  and  triaxial  confinements. 

7.3  VARIATION  OF  STRESS  IN  TWO  PRINCIPAL  DIRECTIONS 

As  shown  in  Fig.  7.1(b),  the  second  series  of  triaxial  con¬ 
finement  tests  started  with  an  Isotropic  condition  of  28  psi  (193  kPa). 
Then,  the  stress  in  the  y-di  recti  on  was  decreased  from  28  psi  (193  kPa) 
to  16  psi  (110  kPa)  in  decrements  of  4  psi  (28  kPa)  while  the  stress  in 
the  z-direction  was  increased  from  28  psi  (193  kPa)  to  40  spi  (276  kPa) 
in  increments  of  4  psi  (28  kPa).  As  a  result,  the  stress  in  the 
x-direction  represented  the  intermediate  principal  stress  which 
remained  constant  at  28  psi  (193  kPa)  while  the  minor  and  major  princi¬ 
pal  stresses  in  the  y-  and  z-  directions,  respectively,  varied  accord¬ 
ing  to  the  pattern  described  above. 

7.3.1  Effect  of  Major  and  Minor  Principal  Stresses 

For  the  stress  conditions  shown  in  Fig.  7.1(b),  P-wave 
velocities  In  the  x-dlrectlon  varied  from  1761  ft/s  (537  m/s)  to  1792 
ft/s  (547  m/s)  with  an  average  value  of  1775  ft/s  (541  m/s)  and  stand¬ 
ard  deviation  of  14  ft/s  (4  m/s).  Therefore,  P-wave  velocities  in  the 
x-direction  (intermediate  principal  stress  direction)  could  be  consid- 


ered  as  constant  at  1775  ft/s  (541  m/s)  and  Independent  of  the  varying 
major  and  minor  principal  stresses.  Consequently,  P-wave  velocities 
in  the  x-direction  appear  as  a  single  point  at  5a  equal  to  28  psi  (193 
kPa)  in  Fig.  7.3. 

Compression  wave  velocities  measured  in  the  major  and  minor 
principal  stress  directions  are  plotted  against  confining  stress  in 
the  direction  of  wave  propagation  in  Fig.  7.3.  The  two  resulting 
best-fit  straight  lines  were  shown  as  solid  lines  in  the  figure.  The 
best-fit  straight  lines  obtained  under  Isotropic  conditions  are  also 
shown  In  Fig.  7.3  as  dashed  lines  for  comparison  purposes.  Since  the 
solid  lines  are  very  close  to  the  dashed  lines,  Eq.  6.1  could  be  used 
to  express  the  relationship  between  P-wave  velocity  and  confining 
stress  In  the  direction  of  wave  propagation  under  this  type  of  trlaxlal 
confinement. 

7.4  VARIATION  OF  STRESS  IN  ALL  PRINCIPAL  DIRECTIONS 

As  shown  In  Fig.  7.1(c),  the  third  series  of  trlaxlal  con¬ 
finement  tests  also  started  with  an  Isotropic  condition  of  28  psi  (193 
kPa).  Then  the  stress  In  the  x-dlrectlon  was  Increased  from  28  psi  (193 
kPa)  to  37  psi  (255  kPa)  In  increments  of  3  psi  (21  kPa)  while  stress 
variations  in  the  y-  and  z-directions  were  the  same  as  in  the  second 


series  of  tests  (Section  7.3).  Therefore,  confining  stresses  In  all 
directions  varied  with  the  stresses  In  the  z-,  x-  and  y-dlrections 
being  the  major,  Intermediate  and  minor  principal  stresses,  respec¬ 
tively. 

7.4.1  Effect  of  Principal  Stresses 


Compression  wave  velocities  obtained  from  stress  conditions 
shown  In  Fig.  7.1(c)  are  plotted  against  the  confining  stress  In  the 
direction  of  wave  propagation  In  Fig.  7.4.  In  this  case,  each  confin¬ 
ing  stress  Is  one  of  the  principal  stresses.  As  shown  in  Fig.  7.4, 
three  best-fit  straight  lines  (solid  lines  in  the  figure)  were 
obtained:  one  for  each  principal  direction.  Best-fit  straight  lines 
obtained  under  Isotropic  conditions  are  also  shown  In  Fig.  7.4  as 
dashed  lines.  Once  again  the  behavior  under  trlaxlal  loading  Is  very 
close  to  that  determined  under  Isotropic  and  biaxial  conditions,  and 
hence  Eq.  6.1  can  be  used  to  express  the  relationship  between  P-wave 
velocity  and  confining  stress  In  the  direction  of  wave  propagation. 


It  should  be  noted  that  due  to  structural  anisotropy,  P-wave 
velocities  In  the  z-dlrectlon  are  less  than  those  In  the  x-  and  y- 
dlrectlons  which  were  approximately  equal.  This  Is  In  same  relative 
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Fig.  7.4  -  Variation  of  Compression  Wave  Velocities 
Under  Triaxial  Confinement  When  All  Prin¬ 
cipal  Stresses  Varied 


ordering  of  wave  velocities  with  directions  as  observed  under  Isotrop 
1c  loading  (Fig.  4.2). 


An  effort  was  made  to  develop  a  generalized  equation  relat¬ 
ing  P-wave  velocities  In  one  principal  direction  to  confining  stresses 
In  all  three  principal  directions  under  trlaxlal  confinement.  This 
generalized  equation  had  the  form  of: 


„  -  r  -  mx  -  my  _  mz 

v  =  C  a  a  v  a 

xx  xx  yy  zz 


(7.1) 


In  which  the  four  constants  had  to  be  determined  from  velocity  and 
stress  data.  It  just  so  happened  that  there  were  four  sets  of  data  to 
fulfill  the  minimum  requirement  of  measured  values  of  variables.  Ide¬ 
ally,  the  power  for  the  stress  In  the  direction  of  wave  propagation 
should  have  a  value  In  the  range  of  0.20  to  0.25  while  stresses  In  oth¬ 
er  two  perpendicular  directions  should  have  powers  close  to  zero.  How¬ 
ever,  results  from  multiple  regression  analysis  had  values  of  the 
parameters  C,  mx,  my  and  mz  far  different  from  the  expected  values  and 
valred  randomly  from  analysis  for  P-wave  velocities  In  each  of  the 
principal  directions.  Therefore,  It  was  believed  that  the  sample  popu¬ 
lation  was  too  small,  and  results  from  multiple  regression  analysis 
were  too  sensitive  to  provide  any  valid  statistical  model  for  the 
P-wave  velocities  and  stresses  In  the  principal  directions  under  trl 


axial  confinement. 
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7.5  SUMMARY 


Three  series  of  trl axial  confinement  tests  were  performed  to 
study  the  effect  of  stress  states  on  P-wave  velocities.  The  results 
from  the  first,  second  and  third  series  of  tests  are  presented  In  Sec¬ 
tions  7.2,  7.3  and  7.4,  respectively.  All  results  leaded  to  the  same 
conclusion  that  It  Is  essentially  the  confining  stress  In  the  direction 
of  wave  propagation  that  controls  the  P-wave  velocity  and  the  relation¬ 
ship  can  be  expressed  by  Eq.  6.1. 

Best-fit  straight  lines  obtained  from  the  average  of  the 
three  series  of  trlaxlal  confinement  tests  are  shown  as  solid  lines  In 
Fig.  7.5.  Best-fit  straight  lines  obtained  under  Isotropic  conditions 
are  also  shown  In  Fig.  7.5  (as  dashed  lines).  Once  again,  both  sets  of 
straight  lines  are  very  close  to  each  other  Indicating  that  Eq.  6.1  can 
be  applied  for  both  Isotropic  and  trlaxlal  confinements.  Constants  for 
these  lines  are  summarized  In  Table  7.1. 


triaxial 

isotropic 


< 


TABLE  7.1 


Comparison  of  Constants  and  Slopes  for  Equations 
Relating  Vp  to  aQ  for  Isotropic  Confinement* 
and  Vp  to  5a  for  Triaxial  Confinement** 


DIRECTION 

Cl 

ISOTROPIC  ,  TRIAXIAL 

m 

ISOTROPIC  .  TRIAXIAL 

X  (NS) 

308 

287 

0.22 

0.22 

Y  (EW) 

280 

282 

0.23 

0.23 

Z  (TB) 

371 

272 

0.17 

0.21 

(Isotropic  Confinement  Eq.  4.2) 


(Triaxial  Confinement,  same  as  Eq.  6.1 
for  Biaxial  Confinement) 


CHAPTER  EIGHT 


SUMMARY  AND  CONCLUSIONS 


8.1  SUMMARY 

One  sample  of  dry  sand  was  tested  In  this  study.  The  sand 
was  the  same  as  that  used  by  Kopperman  et  al  (1982)  and  Knox  (1982)  In 
similar  earlier  studies.  The  sample  was  constructed  with  a  special 
raining  device  (Section  2.3)  constructed  as  part  of  the  project.  With 
this  raining  device.  It  was  possible  to  construct  a  uniform  sample  with 
a  void  ratio  of  0.64  ±0.03  and  dry  density  of  101.8  pcf  ±2.0  pcf  (1629 
kg/m3  ±32  kg/m3). 

Tests  with  various  combinations  of  principal  stresses  were 
performed  with  the  sand  sample,  namely:  isotropic  confinement 
(5i  =  a2  =  03),  biaxial  confinement  (5  >  52  =  o3  or  3i  =  o2  >  o3),  and 
triaxial  confinement  (3j  >  d2  >  53).  In  each  test,  velocities  of  P- 
waves  propagating  along  each  principal  axis  were  measured,  and  the 
results  were  used  to  study  the  effects  of  confinement  on  P-wave 
velocity. 


8.2  CONCLUSIONS 
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From  the  results  of  P-wave  velocity  measured  under  various 
isotropic,  biaxial  and  triaxial  confinement  states,  the  following  con¬ 
clusions  have  been  reached. 

1.  Under  isotropic  confinement,  the  relationship  between 
P-wave  velocity  and  isotropic  confining  pressure  can  be  expressed  in 
the  form  used  by  Hardin  and  Richart  (1963): 

Vp  =  (8.1) 

where:  Vp  =  p-wave  velocity,  in  fps, 

=  Isotropic  confining  pressure,  in  psf, 

Cx  =  constant  (Cx  =  Vp  when  5q  =  1  psf),  and 
m  =  slope  of  the  curve. 

2.  The  effect  of  stress  history  on  P-wave  velocity  is  small 
and  can  be  neglected. 

3.  Due  to  structural  anisotropy,  the  sand  sample  under  iso¬ 
tropic  confinement  can  be  treated  as  a  cross-anisotropic  material  with 
the  horizontal  plane  as  the  plane  of  symmetry.  It  takes  five  independ¬ 
ent  constants  to  determine  completely  the  cross-anisotropic  model. 


Under  biaxial  confinement  with  different  horizontal  confining  stresses 
or  under  trlaxlal  confinement,  stress  anisotropy  couples  with  struc¬ 
tural  anisotropy  to  produce  a  completely  anisotropic  model  for  the  sand 
sample. 

4.  Under  biaxial  confinement,  the  relationship  between 
P-wave  velocity  and  confining  stress  In  the  direction  of  wave  propa¬ 
gation,  5a,  can  be  expressed  by: 

Vp  =  Cxa*  (8.2) 

It  should  be  noted  that  Eq.  8.2  (Instead  of  Eq.  8.1)  can  also  be  used 
for  isotropic  confinement  because  the  mean  effective  principal  stress, 
oq.  Is  equal  to  the  principal  effective  stress  In  the  direction  of  wave 
propagation,  .  Besides,  Eq.  8.2  had  been  verified  for  biaxial  condi¬ 
tions  with  one  or  two  effective  principal  stresses  varying. 

5.  Under  trlaxlal  confinement,  the  same  equation  (Eq.  8.2) 
can  be  used  to  express  the  relationship  between  P-wave  velocity  and 
confining  stress  In  the  direction  of  wave  propagation. 

6.  For  this  sand  sample,  the  general  equation  relating 
P-wave  velocity  and  confining  stress  In  the  direction  of  wave  propa¬ 
gation  can  be  expressed  by  Eq.  8.2,  Independent  of  stress  state.  The 


v.V>, 


Amplitude  Change 


Percent  Strain  Per  Unit  Amplitude  Change 


Percent  Striln  Per  Unit  Amplitude  Change 
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